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SYNOPSIS 



SYNOPSIS 


This thesis presents a study of the systematics of molecular 
motion in liquid systems using nuclear spin relaxation technique. 
Measurements of proton spin-lattice relaxation times, T^ , and 
coefficients of self-diffusion, D, have been made as functions of 
temperature for several liquids from room temperature upto about 
4o°K below the critical temperature or the temperature at which 
they might decompose. A spin— echo spectrometer operating at 
9.78 MHz was fabricated for this purpose. The liquids studied are 
thiophene, furan, 2- and 3-bromothiophenes and 1, 2-dichloroethane. 
The experimental results have been analysed in terms of three 
dominant mechanisms of relaxation, namely, the inter— and intra- 
molecular dipolar interactions and spin-rotation interactions. The 
intermolecular contribution has been separated by using the D 
values measured and the correlation times for molecular reorienta- 
tion have been determined as functions of temperature. In certain 
cases where spin-rotation interaction plays a significant role it 
has been possible to determine the angular velocity correlation 
times, their temperature dependence and estimate the spin-rotation 
interaction constant for the molecules. 

The first chapter of the thesis introduces the concept of 
nuclear spin relaxation with special emphasis on those features 
of relaxation processes in liquids which make this type of study 
a probe of certain characteristics of molecular motion. This is 
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followed by a review of the previous work done in this direction 
and later a discussion of the scope, method and motivation of the 
present work. 

The theoretical framework used in the analysis of the experi- 
mental results of this work is described in Chapter II. The 
equations of motion of magnetization are developed on the basis 
of the density matrix theory to obtain rigorous expressions for 
Tj^ in terms of the correlation functions for molecular motion. The 
theoretical features of various important mechanisms of relaxation 
in liquids have been given. It turns out that only inter- and 
intramolecular dipolar interactions and spin-rotation interactions 
are of interest for the present analysis. A brief discussion of 
Torrey's theory of translational diffusion and Huntress' theory on 
anisotropic rotational diffusion used in the analysis are given. 
The existence of multiple exponentials in relaxation decays in 
molecules containing several spins located at physically inequi- 
valent sites has also been discussed. 

Chapter III deals with the description of the spin-echo 
apparatus fabricated for this work. The essential features of 
this apparatus are a pulse programmer comprising of Tektronix 
waveform and pulse generators and some simple circuits. This is 
used to program the pulses in a desired sequ.ence. These d.c. 
pulses are used to gate a continuous-wave crystal-controlled 
oscillator operating at 9,78 MHz. At the output of the gated 
oscillator radio-frequency (r.f.) voltage of about 300 volts peak 
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to peak is obtained. The r.f. receiver is a five-stage amplifier- 

detector assembly tunable from 8-2o MHz and has a net gain of 

about 10,000 and a bandwidth of 1 MHz. It has a signal to noise 

ratio of 2:1 for 8o% modulation at 3 [i v carrier voltage and takes 

a maximum input signal of 300 M’ v. Another important feature of 

the apparatus is a variable temoerature cryostat designed to 

measure relaxation times from liquid nitrogen temperatures up to 

about 6oo‘^K. All the experiments in this work are performed along 

the liquid-vapor equilibrium curve by appropriately filling the 

sample tubes. The samples are prepared free of dissolved oxygen 

by using freeze-pump -thaw technique and gettering. 

Chapter IV describes a study on thiophene and furan in the 

temperature range 300°K to 550°K and 295°K to 4o5°K respectively. 

Tj_ varies from 24.0 sec at 300 °K to 7 2,0 sec at 55o°K with a 

maximum of 94.0 sec at 490 °K for thiophene. The coefficient of 

—5 2 ««5 

self-diffusion varies from 2,4 x lo” cm /sec to 38.0 x lo” 

cm /sec in the same temperature range. In the case of furan T^ 

varies from 39.0 sec at 295°K to 66.0 sec at 4o5'^K with a maximum 

of 95.0 sec at 373‘^K and D varies from 5.5 x lo"*^ cm^/sec to 
-5 p 

36.0 X 10 cmvsec. The contribution of the intermolecular 
dipolar interaction has been separated using Torrey' s theory. 
Analysis of the intramolecular interactions gives the temperature 
dependence of , the rotational diffusion constant perpendicular 
to the plane of the molecule and Tj , the angular velocity 
correlation time. Di varies from 0,21 x lo sec at 300°K to 



xii 


3.5 X 10^^ sec“^ at 550°K in thiophene and from 0.37 x sec ^ 

at 295°K to 4.02 x lo^^ sec“^ at 4o5‘^K in furan. r j varies from 
1.4 X sec at 300^K to 10.8 x lo”^^ sec at 550°K in thiophene 

and from 1.6 x lo"^'^ sec at 295'^K to 12.7 x lo"*^"^ sec at 4o5°K in 
furan. Spin-rotation interaction constants for thiophene and 
furan have been estimated to be about l.o kHz and 1.4 kHz respect- 
ively. The reorientation is found to be diffusional in character 
upto 48o°K in thiophene and 380 °K in furan. The variation of 
with temperature in thiophene obeys an Arrhenius type of equation 
upto 500°K with an activation energy of 3.2 k cal/mole. However, 
it does not show such a behaviour in furan. 

Chapter V presents T^ and D measurements on 2- and 3-bromo- 
thiophenes in the temperature range 300 °K to 45o°K and their 
analysis. Higher temperatures were avoided due to incipient 
decomposition of the samples. The results show indication of 
only dipolar interactions contributing to relaxation. The analysis 
of these results, though somewhat crude due to lack of data on the 
structural and motional constants, reveals an interesting contrast 
in the temperature dependences of the rotational correlation times 

for these molecules. y~ varies from 0.23 x 10~^^ sec at 333®K 
11 o 

to 0.11 X lo~ sec at 450°K with an activation energy of 2.0 

k cal/mole in 2-bromothiophene while^in 3 -bromo thiophene it remains 

—1 1 

essentially constant at about 0.7 x 10 sec over the entire 
temperature range indicating a large activation energy. 
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Chapter VI deals with a similar study on 2-dichloroethane 
in the temperature range 300°K to 430°K. The possibility of 
scalar coupling of protons with chlorine nuclei and intramolecular 
dipolar interactions modulated by rotamer exchange contributing 
to relaxation has been examined and found to be negligible. The 
rotational correlation time is found to vary rather slowly with 

p 

temperature. This variation is similar to the T^ dependence 
predicted by Steele for a spherical top undergoing reorientation 
dominated by inertial effects. Confirmation of this interpreta- 
tion^ however, needs further experimental and theoretical work 
on similar molecules. 

The significant results of this work are summarized in 
Chapter VII. 

The thesis has three apoendices: (i) a note on the reciprocal 
temperature dependences of the coefficient of self-diffusion, D, 
and coefficient of viscosity, r) , in several organic liquids, 

(ii) a collection of some of the simple circuits used in the 
apoaratus which are not included in Chapter III and (iii) a 
description of the Carr-Purcell method for T 2 measurement and 
its adaptation to the spectrometer. 
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INTRODUCTION 


I .A Nuclear Spin Relaxa-tion in Liquid Systems 


Consider an assembly of nuclear spins each of 
angular momentum I In the absence of any external magnetic 

field the energy levels are (21+1) -fold degenerate and they 
are, therefore, equally populated. Now, if a static magnetic 

A 

field Hq^c is applied this degeneracy is lifted to give a set 
of (21+1) equidistant levels with energies (-ynliHo) where 
m = -I, -I+l, .... +1. When the field is switched on the 
system is in a non-equilibrium condition in which these levels 
are equally populated. The spins eventually attain thermal 
equilibrium represented by a Boltzmann distribution of popula- 
tions (at the temperature of the surroundings) among the energy 
levels by dissipating the excess energy to the other degrees of 
freedo-no-. The process by which this equilibrium is attained is 
bnown as relaxation. This is the same process by which a 
magnetization = Xq Hq (Xq is the static susceptibility) is 
obtained along the direction of the field and any magnetization 
in the perpendicular plane is made to vanish. The motion of 
the magnetization M is represented by a phenomenological equa- 
tion, given by Bloch^, 


dM 

dt 


= YMx H 



1c 


Mg - Mq 
Tl 


(I.l) 


where , My and Mg represent the instantaneous magnetizations 
along the x,y and z axis respectively. The first term refers to 
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the precession of the magnetization vector ajoout the z— axis 
with the Larmor frequency 

Wq =YHq * (1.2) 

governs the growth of magnetization along the z-axis and is 
called the spin-lattice (or longitudinal) relaxation time and 
T2 the decay of magnetization in the xy-plane and is referred 
to as the spin-spin (or transverse) relaxation time. These 
two times can, in general, be different since the growth of 
populations of the spin states (given by T^) can proceed at a 

rate different from that at which the phase coherence between 

, 2 
these states is destroyed (given by T2) . 

In a nuclear magnetic resonance (nmr) experiment a radio- 

frequency (r.f.) magnetic field at the Larmor frequency induces 

transitions between the energy levels causing a net absorption of 

3 4 

energy by the spin system ' . The strength of the r.f. field 
that can be used in order that a sustained absorption signal is 
obtained is governed by the efficiency of relaxation. 

Relaxation processes arise due to interactions within the 
spin system (spin-spin interactions) or due to interactions of 
the spins with other degrees of freedom termed as 'lattice* 
(spin-lattice interactions) . These interactions manifest 
themselves in different ways depending upon the physical state 
of the sample. In diamagnetic solids at low temperatures, for 
example, the dipole-dipole interactions between nuclear spins 
produce local fields of the order of several gauss resulting in 
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broad resonance lines and allows the spin system to attain 
internal thermal equilibrium much quicker than with the lattice 
i.e, the T^s are much shorter than T^s. In a description of 
relaxation in solids it is necessary to consider all the nuclear 
spins as forming the spin system because of the tight coupling 
between them"^ ' ^ . 

In the case of liquids (and gases) / however, the situation 
is quite different. Due to the rapid molecular motion the 
local fields are usually averaged out leading to narrow reso- 
nance lines (having linewidths of the order of fraction of a 
milligauss) for spin-Jg nuclei in liquids^. This also results 
in T^s of the same order of magnitude and often equal to T^s 
in these systems. In the description of relaxation in lic[uids 
it is sufficient to consider a group of spins on the molecule 
(sometimes even the individual spins) as the spin system"^. The 
molecular motion plays a crucial role in determining the effi- 
ciency of spin relaxation in liquids and, therefore, the study 
of relaxation under chosen physical conditions could be a 
significant means for the study of some aspects of molecular 
motion , In the next section this point is discussed in some 
detail. 

It may often happen that the relaxation of the spin system 
does not follow the simple pattern set forth by Eq. (I.l). For 
solids Eq, (I.l) is found to be certainly inadequate . In the 
case of liquids, however, the form of Bloch equations is found 
to obtain under a number of generally readiatable conditions 
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although there are exceptions, A rigorous description of the 
relaxation effects may be obtained by using the density matrix 
method and this theory indicates clearly the cases where 
deviations may be expected from Eg. (I.l) ' . The features of 

this theory, in so far as they are pertinent for the analysis 
of the present work, are summarized in Sec. II. A, 


I.B Spin Relaxation as a Probe for Molecular Motion 


The understanding of nuclear spin relaxation in 

liquids is chiefly due to the classic theory of Bloembergen, 

14 

Purcell and Pound (BPP) . It was mentioned earlier that in a 
liquid the local fields at the site of any nucleus, arising from 
other nuclei in the sample, are averaged out because of the rapid 
molecular motion. In this theory it was shown that the effi- 
ciency of spin relaxation is given by the spectral density, at 
the nuclear Larmor frequency, of the time correlaticsai functions 
of such rapidly fluctuating fields. Any spin-lattice interaction 
may be written 

H'(t) = \ (-1)^ F‘5(t)A“^ , (1,3) 

q 

where P^(t) are the lattice variables which are random functions 

of time and A^ are spin operators. The relaxation rate due to 

7 14 

this interaction is determined by ' 

J(a>) = /“ e-'^“'^ G(t) dm, ( 1 . 4 ) 


where 


G(t) = I G (t) 
q ^ 


I<F^*(t)F^(t+1^ 

q 


av 


l 

q 


<F'3* (o)Fq(T;)>^^ ^ 

(1,5) 
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with < > _ representing an ensemble average. F'^(t) are 

considered stationary random functions of time, G(t: ) is the 
time correlation function and J(w) is its spectral density at 
the angular frequency oj. If F'^Ct) is given by a Markoff process 
the correlation function decays exponentially: 

< F^*(o)F^(T^) > =<1 f'^( 0) l^>e“ , (1.6) 


Where is the correlation time for the appropriate molecular 
motion. Substituting Eq. (1,6) in Eq. (1.4) gives 


Tl q 


I <jF^(0) 1 > 




4v/ 1+ 


(1.7) 


It may seem from Eq, (1.4) that a determination of the 
relaxation rate, arising from an interaction of known strength, 
as a function of u would allow, by Fourier inversion, the deter- 
mination of G(t ) as a function of t for the appropriate 
lattice function. However, this is not normally feasible since 
the Larmor frequencies used in nmr experiments range between 
1 to loo MHz while the motional correlation times are in the 
range of lo“^^ to lo"^^ sec so that w « 1. This results in 
the spectral density being independent of frequency, a condition 
known as 'extreme narrowing'^. The relaxation rates, therefore, 
provide information on the correlation time given by the 
integral of the normalized correlation function 

r<pq*(0)F^(T) > dt 


'"c = 


( 1 . 8 ) 


'<w 


It can be seen that the relaxation times are not sensitive to 
the actual "shape" of the correlation function but depend on the 
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areas enclosed by these functions. Such integrals of correla- 
tion functions can also be related, quite generally, to appro- 
priate macroscopic transport coefficients^"^ 

We shall now consider the various types of correlation 
times that arise in the description of relaxation processes in 
liquids. It will be seen that for proton spin relaxation in 
diamagnetic liquids studied in this work the only likely mecha- 
nisms which contribute significantly to relaxation are the inter- 
and intramolecular dipolar interactions and the spin-rotation 
interactions. Why these are the only ones shall be discussed 
in some detail in Sec. II. B. The intermolecular dipolar interac- 
tion is primarily modulated by the relative translational motion 
of the molecules which involves a translational correlation time 
T-t- The intramolecular dipolar interaction is, on the other hand, 
modulated by the reorlentational motion of the molecules involv- 
ing a rotational correlation time The spin-rotation interac- 

tion, which is a magnetic interaction of the nuclear spin with 
the rotational angular momentum of the molecule * , is modulated 

by changes in the angular velocity of the molecule as it 

O O 

reorients . The corresponding correlation time is the angular 

velocity correlation time t j, in certain favourable circumstances 
it is possible to separate the contribution coming from the 
different mehcanisms, whence the correlation times of molecular 
motion t -j- and Tj can be estimated. The theoretical calcula- 

tion of the correlation functions for these interactions is 
discussed in Sec. II. C and it is shown that if the translational 
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and reorientational motions are diffusional in character t ^ is 
related to the coefficient of molecular self-diffusion and-r ^ 
and to the elements of the rotational diffusion tensor. 

The next section summarizes some of the previous theoretical 
and experimental work aimed at extracting such information on 
molecular motion in liquids through spin relaxation and in 
Sec. I.D the scope of the present work is discussed. 

I.C Previous Work on Molecular Motion by Spin Relaxation 

14 

With the BPP theory providing the basic formulation 

of the relaxation rates in terms of correlation functions and the 

13 

density matrix theory of Wangsness and Bloch the rigorous form 
of the differential equations for macroscopic magnetization, the 
theoretical problem of spin relaxation in liquids becomes that 
of calculation of the appropriate correlation functions (or 
correlation times) on the basis of specific descriptions of 
molecular motion in liquids. 

In the original BPP theory the correlation time-t^ for 
intermolecular dipolar interaction modulated by translational 
diffusion is written in terms of the coefficient of self -diffusion 
D and a hard-sphere radius 'a* for the molecule = a^/12D) , 

and D was related to the coefficient of viscosity T] through 
the Stokes -Einstein relation (D = kT/6tT)a, k is Boltzmann 
constant and T is the temperature in ®K) . Torrey^”^ has 
considered this problem in a somewhat rigorous fashion on the 
basis of the theoiry of random flights which, in the limiting case 
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of small mean squared flight paths and small times between 

flights, goes over to the results given by the equation for 

isotropic translational diffusion. The derivation of the 

correlation time in terms of the diffusion equation is briefly 

discussed in Sec. II, C. For the case of spherical molecules, 

each containing a spin located at its center, this correlation 

28 

time is the same as that given by BPP. Hubbard considered 

the case of molecules in which the nuclei are not at the center 

and showed that this leads to a small correction to the above 

result. Torrey's formula leads to fairly correct estimates of 

relaxation rates arising from intermolecular interactions. 

For the reorientational motion BPP have used the Debye 

model of isotropic rotational diffusion originally used in the 

29 

theory of dielectric relaxation in polar liquids . This model 
assumes that the medium surrounding the molecule is continuous 
so that the rotational diffusion constant D^. is given by the 
Stake's fo3rmula kT/Sn T) a. There is considerable evidence that 
the contribution to relaxation rates calculated by \ising this 
formula is much larger than that observed experimentally^^' 

The major shortcomings of Debye model are (1) the rotational 
diffusion may, in general, be anisotropic (2) there may be 
significant contribution of inertial effects to reorientational 
processes (3) the use of Stokes -Einstein relation between Dj- and 
11 may be questionable for most liquids except possibly for large 
molecules in solution. 
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The problem of anisotropic rotational diffusion has been 

treated by Favro^^z Huntress^^z Ivanov^^ and others^^”^^. In 

1 R 

the analysis of the results of the present work Huntress* the'ory 

has been used and shall be discussed in some detail in Sec. II. C. 

This theory was earlier used to considerable advantage in the 

analysis of quadrupolar relaxation data to determine the princi~ 

pie values of the rotational diffusion tensor in several 
37—40 41 

liquids . Steele has evaluated cdrrelation functions where 

inertial effects are predominant in molecular reorientation 
starting from the rotational analogue of Langevin equation. 

Some aspects of this are given in Chapter VI. The approximation 
of the medium as a continuum, implicit in Stokes formula, is 
discarded in the above. 

The spin-rotation interaction as a mechanism for relaxation 

23 24 

in liquids has been discovered recently ' . In gases, however, 

this mechanism has been known to contribute significantly to 
• 42 

relaxation . Here the molecular collisions that change the 

orientation also change angular velocity (or rotational angular 

momentxam) and therefore and are the same (in the strong 

collision limit). In the case of liquids, on the other hand, 

the two correlation times are quite different since angular 

velocity occurs only during a molecular collision and the 

correlation functions are different from those for orientation 
24 

changes . The existence of this interaction for fluorine 
relaxation has been observed'^^ ' and recently it has also been 
seen to contribute to proton relaxation at high temperatures by 
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Powles and coworkers . It is important to note that 

the tempera tute dependence o£ the relaxation from spin-rotation 
interaction is opposite to that from dipolar interaction and 
hence the two could be separated if the temperature dependence 
of relaxation is determined. Another advantage of such a study 
is that it enables the determination of spin-rotation interaction 

constants of small magnitudes ( 'v llcHz ) which is not usually 

45 

possible by other methods like molecular beams and microwave 
21 

spectroscopy . The correlation times for this interaction can 

be obtained by using the theoiry of anisotropic rotational 

15 

diffusion given by Huntress 

There is at present, no general theory available for 

reorientational motion of molecules which can be applied to all 

situations. While the limiting conditions of diffusional motion 

assuming small angular steps in the rotational random walk and 

that of inertial motion implying essentially free rotation can 

be handled the intermediate situations are intractable. 

Furthermore, in any description of reorientational motion the 

molecular shape plays an important role and it is difficult to 

simplify the theory to a usable form for all molecular structures. 

Considerable experimental evidence exists from neutron scattering 

46 

measurements for the need of a general theory , and the results 
of nuclear spin relaxation, though somewhat limited in scope by 
comparison ^ throw considerable light on the need for competant 
understanding of the reorientational processes in liquids. 
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From the experimental point of view the important problem 
is to separate the individual contributions of the different 
interactions responsible for relaxation in a liquid. An experi- 
mental method for explicit separation of the inter- and intra- 
molecular interaction for proton relaxation is to study the 
variation of relaxation time as a function of concentration in 
mixtures of the liquid and its deuterated form'^ * * . Such 

separations have also been affected by using experimental values 

of T 2 _ and known formulae for inter- and intramolecular interac- 
26 

tions . In the present work and D values measured as 
functions of temperature are 'analysed to affect such a separation. 
The spirit of this work is described briefly in the next section. 

I.D Scope of the Present Work 

The work presented in this thesis is based on an 
experimental study of spin-lattice relaxation as a function of 
temperature in five different liquid systems (i) thiophene 
(ii) furan (iii) 2-bromothiophene (iv) 3-bromothiophene and 
(v) 1, 2-dichloroethane and is aimed at obtaining the character- 
istics of molecular motion in these liquids. The experiments 
were performed on a spin-echo spectrometer operating at 9,78 MHz 
fabricated for this purpose. The temperature for different 
liquids ranged from room temperature upto -v 40 °K below the 
critical temperature or the temperature at which they might 
decompose. 
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In the analysis of the data the contribution of the inter- 
molecular dipolar interaction is separated from the observed 
relaxation rate by explicitly determining the coefficient of 
self-diffusion in situ ^of' the sample as a function of tempera- 
ture. An accurate measurement of D is afforded easily by the 

4-9 

spin-echo method’ . The use of some additional data on density 
and viscosity along with D values allows a fairly correct 
estimate of this contribution. 

All the molecules studied contain several spins all of 
which do not necessarily experience intramolecular dipolar 
Interactions of the same strength. When this feature is included 
in the general theory of relaxation by dipolar interaction one 
finds that the T^-decay is usually given by multiple exponentials 
due to cross-relaxation effects * . The actual decays observed 

in the experiments are, however, single exponentials. This is 
explained by the fact that the other time constants have ampli- 
tudes too small to be observed. Nevertheless, a detailed consi- 
deration of the cross -relaxation effects in the present analysis 
led to expressions for relaxation rates that are not linear 
superpositions of the contributions of the individual mechanisms. 

The intramolecular contribution to relaxation obtained in 
this manner is analysed to determine the reorientational para- 
meters in these molecules as a function of temperature. The 
salient features of these results are given below. 

The experimental results on thiophene and furan fall in 
a similar pattern and show a maximum in the T2_ vs lo /T curve 
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Which is a clear indication of spin-rotation interaction being 
present as a mechanism of relaxation for both these compounds. 

The analysis of the intramolecular relaxation rate in terms of 
the theory of anisotropic rotational diffusion is simplified 
by assuming that Tj , the angular velocity correlation time, 
is isotropic. This relates the three principal diffusion 
constants through the respective moments of inertia. With this 
assumption and the known data on molecular structure, the 
temperature dependence of the rotational diffusion constant 
perpendicular to the plane of the molecule (both these molecules 
are planar) and Tj are determined for both molecules. Furthermore, 
the spin-rotation interaction constants have been estimated. 

The results indicate that the molecular reorientation in these 
systems show deviations from diffusional character at high 
temperatures . 

3 

For 2- and 3-bromothiophenes the Tj_ vs lo /T curves, however, 
give no indication, of the presence of spin-rotation interaction. 
The intramolecular relaxation rate is, therefore, entirely due 
to internal dipolar interactions. The temperature dependence of 
the rotational correlation times has been derived for both these 
molecules and are seen to be quite different from each other. 

The analysis of the data on these molecules, however, required 
several assumptions regarding the structural and motional const- 
ants since these are not available in literature. The contrast 

in the variation of for the two cases is, even after allowing 

(s 

for this, unmistakable. 
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In 1, 2-diGhloroethane the intramolecular relaxation rate 
shows a very small variation with temperature. The reorienta- 
tional motion of this molecule has the additional complication 
that the CH 2 CI groups undergo hindered rotation about the C-C 
bond while the entire molecule tumbles. The rotamer exchange, 
however, does not seem to be an effective mechanism for modula- 
ting the dipolar interactions. The variation of relaxation 
rates with temperature indicates that the molecular reorienta- 
tion in this case is predominantly inertial in character. 

The material in this thesis is organized in the following 
mannert The theoretical frameworlc for the analysis of experi- 
mental results is sketched in Chapter II. This includes the 
development of macroscopic differential equations on the basis 
of the density matrix theory, features of the relaxation 
mechanisms in liquids, derivation of the relevant correlation 
functions of molecular motion and a discussion of cross- 
relaxation effects leading to multiple-exponential relaxation 
decays. 

The spin-echo apparatus fabricated by the author for this 
work is described in Chapter III. This contains the principle of 
the method, discussion of the essential circuitry, performance 
of the equipment, methods of measurement of relaxation times and 
diffusion constants, a description of the temperature variation 
accessory and details of sample preparation. 

Chapter IV, V and VI contain the experimental results and their 
analysis, on thiophene and furan, 2- and 3 -bromo thiophenes and 
1, 2-dichloroethane respectively. 
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The important conclusions of this work are summarized 
briefly in Chapter VII. 

This thesis has three appendices/ namely/ (i) Temperature 
dependence of the coefficient of self-diffusion in liquids/ 

(ii) Additional details of the spin-echo apparatus/ and 

(iii) Measurement of T 2 by Carr -Puree 11 method. 



17 


Chapter II 


THEORETICAL FRAMEWORK 

Page 

II. A Density Matrix Theory; Macroscopic 

Differential Equations 18 

II. B Relaxation Mechanisms in Liquids 23 

(1) Dipolar interaction 23 

(2) Spin-rotation interaction 24 

(3) Quadrupolar interaction 26 

(4) Scalar coupling modulated by 

(a) Chemical exchange 26 

(b) Rapid relaxation of spin S 26 

mm 

(5) Anisotropic chemical shift 27 

II. C On Correlation Functions 27 

(1) Translational diffusion; 28 

Intermolecular dipolar interaction 

(2) Rotational diffusion; 32 

Intramolecular interactions 

(a) Intramolecular dipolar interaction 33 

(b) Spin-rotation interaction 35 


II. D 


Equations of Motion for Multispin Systems 


38 



18 


THEORETICAL FRAMEWORK 

In this chapter the basic theoretical framework used in 
the analysis of the experimental results obtained in this work 
is described. In Sec. IIA the equations of motion of the 
magnetization are developed on the basis of the density matrix 
theory to obtain rigorous expressions for T^ in terms of the 
correlation functions of molecio^lar motion. Section II. B sketches 
the theoretical features of the important mechanisms of relaxa- 
tion. The calculation of the relevant correlation functions on 
the basis of specific models for molecular motion, along with 
expressions for the correlation times used in later chapters is 
given in Sec. II. C. Finally, in Sec. II. D the case of multi- 
spin systems in which different groups of spins experience 
intramolecular dipolar interaction of different strengths, lead- 
ing, in principle, to multiple-exponential relaxation decays, is 
considered, 

II. A Density Matrix Theory; Macroscopic Differential Ecruations 

A brief discussion of the density matrix formalism 

I 13 

for nuclear induction due to Wangsness %nd Bloch is given 
below in order to develop the equations of motion of the magnet- 
ization and expressions for the relaxation rates. The procedure 

7 

followed is essentially that given by Abragam and what is given 
below is for the puirpose of introducing the notation and some of 
the relevant physical assumptions. The equation of motion for the 
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density matrix p of the system (spin and lattice) is given by 

If = -i[ «o ^ ^ P 

where P and f( ' refer to the spin hamiltonian, the 

hamiltonian of the lattice and the spin-lattice interactions 
all expressed in angular frequency units. Transforming to an 
interaction representation such that 


H'*(t) 


gives 


jiHo^ e^^^ H* - e^^o^ H'(t)^^o^ 


P (t) 


i K^t „iFt ^ ^-iPt „-i ^ 


P e 


(II. 2) 
(II. 3) 


★ 

d p 

dt 


-i [ tf ' * ( t) , P* ] . 


(II. 4) 


Integrating Eq. (II.4) in two successive approximations and 
differentiating with respect to time we have 



= -i[ H'*(t), 


P*(0)] - 


/[W'*(t) I H’*(t-T) 


p (o)] ] du. 


(II. 5) 


Higher order terms are neglected assuming the relaxation to be 
weak. If the lattice is assumed to be in thermal equilibrium 
at t = o , P (o) is factorizable into spin and lattice parts. 

A spin density matrix cr is defined through the relation 

a = Trf {p} and q * = Trf{p *}, (II. 6) 

where Trf { } represents a trace over the lattice states. 

Equation (II, 6) implies 


i H^t -i 
= e ^ ere 


Hot 


cr* 


(II. 7) 
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This leads to 


dt 


- = - Trj { / '^ a-J [K ’*(t) , [ , X* a*(0)]]). (II. 8) 


★ 

is the equilibrium density matrix for the lattice 


states given by 


X _ exp ( - 1iF /TcT) 

o Trf {exp(- TiF/kT)) " 


(II. 9) 


and is the equilibrium spin density matrix at the lattice 


temperature given by 
a*(0) = 


exp(- liH^/hT) 


(II. 10) 


® Tr [exp(- ti HqAT)]' 
where k is the Boltzmann constant. The trace over the lattice 
states for the first term in Eg, (II. 5) gives zero due to the 
particular choice of the time t = O, Equation (II. 8) can be 
'further simplified by using the following considerations : 

(i) In the limit of short correlation times d (o) can be 

k 

replaced by d (t) , 

(ii) For T ^ « t « T^ or T 2 the upper limit of integration can 
be extended to «. 

(iii) li^ile taking Trf{ } high -temperature approximation is made 
for the spin system which leads to d*(t) being replaced by 

(d (t) -C q) and coming out of the double commutator. 

(iv) For liquids the lattice can be described classically and 

Trf { } after multiplying by X^can be replaced by an ensemble 

average <> . This leads, in the Schrodinger picture, to 

dd 
dt 


= -i[HQ,q]-- /dT<[H’(t),[ ff' (t-T)e^”o^ . d> d^]]> 


iH T 


O' - av 
(II. 11) 
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In liquids containing several spins of the same species 
the spin hamiltonian contains the Zeeman interaction of 
all the spins including their chemical shift's, and indirect spin- 
spin interactions. This leads to fine structure in high 
resolution nmr spectra^. In the present study (using spin-echo 
method) the nuclear spins are perturbed by strong r.f, pulses and 
the results may be adequately described by neglecting the 
chemical shifts and soin-spin coupling constants, such that 


Ho = - Wo I • (11.12) 

It has been stated (see Sec, I.B) that the spin-lattice 
interaction hamiltonian can be written as 


where 


(t) = I (-1)^ P^(t)A"^, 
4 


= (- 1 )"^ , 
(A^)^= (-1)^5 


(11.13) 


(11.14) 


where the spin operators are expressed in their irreducible form. 
From Eqs. (II, 12) to (11.14) one gets 




<3/ ^ q iq w t 
A (t) = A e ^ f 


(11.15) 






= I (-l)'" F'"(t) A 
g 


-q -iq » T 
e ^ , 


With this form for H' (t) the correlation functions simplify 
through the relation 

<F^*(t) P^' (0)> < F'^*(t)F^(0)> 

a V 


av 


(11.16) 
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r-j 

This Is Tcnown as the "secular approximation" . Equation (II. 11) 
then becomes 

H = -i[«c„CT]- _ (11.17) 

q. t 

where 

j'^Cq e ) = /“dT <F^ (t)F'^(t-t)> (II. 18) 

° ° av 

The first term in Eq, (II. 17) represents Larmor precession 
and the second term the effect of relaxation on the evolution 
of the system. For the purpose of calculating the effect of 
relaxation the first term need not be considered. 

For a given spin operator Q the macroscopic observable is 
given by 

< Q > = Tr { cT (t) Q}. (11.19) 

The e(quation of motion of < Q > is (from Eq, (11,17)) 

- I (-d'I j'5(qi.^)Tri[A'3,[A‘TQ]] (a-o-oU <11-20) 

For example/ the equation of motion for <I„> is given by 

= - I (-1)'^ j'5(q»o)Tri[Ai [a"^, 1^1] (q-q^)}. 

( 11 . 21 ) ‘ 

7 q 

In the "extreme narrowing" limit u (q u^) may be replaced 
by J(o) , independent of both q and , leading to 

= - I (-1)'* J (0)Tr{[A® [a"! Il](q_ab)}. 

dt q ^ . 


(II. 22) 
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If the simplification of Eq, (II.22) for a particular 
mechanism gives an equation of the form 

- U(<i^> - i^), (II. 23a) 

where = Tr (l^ 0^) ^ the validity of Bloch equations is 
verified for that particular mechanism with 

~ = AJ(0) , (II. 23b) 

where A is a constant of proportionality . 

II.B Relaxation Mechanisms in Liquids 

In this section the various mechanisms which might 
cause relaxation in liquids are described. It will be seen 
from the discussion that only some of them are relevant for 
the liquids investigated. 

( 1) Dipolar interaction ; For a group of spins the 
hamiltonian for this interaction is given by 

H' (t) = I I(i) -D. , (t) •!( j) , (11.24) 

i< j - -IJ 

where D. . is the dipolar interaction tensor between nuclei 

wX J 

i and j. The lattice functions are given by 

1.1=1 Y. Y. 

=.(4fl/5r^Jl (11.25) 

Where r^^j is the intemuclear vector and (o^j / 9j^j) are its 
polar angles in a space-fixed frame. The spin operators are 
given by 
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AT^ =+tl (l + (i)l 2 (j) + 13 ( 1 ) I+(j)) (11.26) 

A?. = (3I^(i)I„(J) - I(i).I(j)), 

where 

+ ily . (11.27) 

There can be two types of dipolar interactions; (a) Interactions 
between two dipoles sitting on the same molecule or Intramolecular 
dipolar interactions. In this case only the reorientation of the 
molecule need be considered and this modulates the orientation 
( 0 , 9 ) of fho internuclear vector. The correlation time involved 
is the rotational correlation time Cb) Interaction between 

dipoles sitting on two different molecules or intermolecular 
dipolar interactions. Here the translational motion of the 
molecules with respect to each other is of relevance and it is 
the magnitude of the internuclear vector r that is modulated - 
because of molecular motion. The correlation time involved 
is the translational correlation time 

This mechanism of relaxation gives an important contribution 
to T^ in most liquids. 

(2) Spin-rotation interaction : This arises due to the 
interaction of the nuclear spin with the magnetic field produced 
at its site by the rotation of the molecule. It represents a term 
bilinear in the coupling of the'VtfepjiwiBtangular momentum of the 
molecule with the angular momentum of the electrons and the 
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coupling of the angular momentum of the electrons with the 
nuclear soins in a second order perturbation calculation of the 
energy of a rotating molecule. If the molecular rotation is 
considered classically this interaction is modulated by changes 
in the angular velocity of the molecule and the correlation time 
involved is the angular velocity correlation time ' . The 

form of the interaction hamiltonian is 


H’(t) = I.c.j(t) = I lic^«j’(t) = y (- 1 ) a' , ( 11 . 28 ) 




with 


where 


= I (-1)“^ C(1U, m-qqm) C'^ (II.29) 

^ iil ^ 


C'° = -1- Tr C* 

o 3 

=r = v'2 




■yz 




(II. 3o) 


■1 

. ( 1 ) 


- 


'"xi = ^ (c — C ) 
IJ 2 


c;( 2 )d,(c;j + ^ Tr c'Sij) 


ij 


and 


A‘° = l' 


(11,31) 


,•+1 ^ T -' ■ 


In Eg. (II.28) I refers to the nuclear spin^C is the spin 
rotation interaction tensor and J is the angular momentum of 
the molecule. Primed quantities are with respect to a molecule - 
fixed coordinate frame. are components of the irreducible 
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Spherical tensor CJ, of ranks. (=0,1,2) contained in C* C^j are the 
cartesian components of C. C(lljL,m-q qm) is a Clebsh-Gordon 
coefficient. In certain cases this interaction forms an 
important mechanism of relaxation and will be considered later. 

(3) Quadrupolar interaction : The interaction of the 
nuclear quadrupole moment with the fluctuating electric field 
gradients at its site is a powerful mechanism of relaxation of 

7 

nuclei with spin > in liquids . Since the present v7ork is 
confined to proton spin relaxation this interaction need not be 
considered any further. 

(4) Scalar coupling modulated by 

(a) Chemical exchange? This is a bilinear coupling 
of the type AI.S in which the spin S undergoes chemical exchange 
If the chemical exchange time Iq much smaller than 1/A this 

7 

coupling becomes a mode of relaxation . The coupling constant 
Jumps randomly between values A and zero corresponding to I and 
S being on the same or on different molecules respectively. 

There is, however, no reason to believe that such a mechanism 
contributes to relaxation in the liquids under study, 

(b) Rapid relaxation of spin §; A situation may 
also arise when the coupling between the spins AI.S becomes time 
dependent when the spin S has a strong mechanism of relaxation 

m- 

of its own such that its relaxation time is much shorter than 
1/A, This often arises when S is a quadrupolar nucleus. The 
spin I, therefore, sees a field produced by S which is a rapidly 
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fluctuating function of time. This is pertinent in 1/2- 

dichloroe thane due to the coupling between protons and chlorine. 

However, it can easily be shown that under normal experimental 

7 , 

conditions this would contribute only to T 2 and not to T]_. We 
could, therefore, ignore this mechanism also in the analysis of 
data of 1, 2-dichloroethane, 

( 5 ) Anisotropic chemical shift : The Zeeman coupling of 
a spin I with a dc magnetic field H given by - sometimes 

requires a coirrection term because of chemical shielding. This 
is given by - H.A.I where ^ is the chemical shielding tensor. 
Because of the molecular tumbling the anisotropic part of the 
shielding tensor contributes to relaxatioit^. There is, however, 
no unambiguous evidence, so far, that such a mechanism contributes 
significantly to relaxation, 

II. 0 On Correlation Functions 

In this section the theoretical calculation of 
correlation functions used in the analysis is described. Prom 
Eqs. (11,18) and (ll,23t)) it can be seen that the relaxation 
times are given in terms of time correlation functions of 
stationary random functions of the lattice variables. The time 
correlation function of P'^CxCt)) is defined as 

G(t) = <F^ (t)F'^(t-K)>av 

= // (x^, t)P'^(x 2 / t+'r)p(x^, t)P(xj_, t;X 2 , tHi-'r)dx^ '^ 2 * 

( 11 , 32 ) 
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where x is the lattice variable which changes with time, p(x 2 ^) 
is the probability of its possessing a value x^ at time t and 
P(x]_, t;X 2 / t+T) is the conditional probability that if x has a 
value X]_ at time t, it would have a value X 2 at time (t+t) . 
pCx^) is easily obtained by considering a completely random 
distribution for x at the initial time t (GCt) does not depend 
on t) . The problem of calculating G(t) then depends on obtain- 
ing an expression for the conditional probability P(x^/ t;X 2 / t+'c) . 
As mentioned before the correlation functions required in 
the present work involve relative translational motion, reorienta- 
tional motion and angular velocity changes of the molecules. 

All these will be considered below in the appropriate diffusional 
limit. 


(1) Translational diffusion ; 

Intermolecular dipolar Interaction : Torrey has given 

a rather general theory of spin relaxation through translational 

27 

diffusion on the basis of the theory of random flights . In 
the limit of small mean flight paths the results of this theory 
approach those given by the translational diffusion equation 

“T* P(r,t) = Dv2p(r,t) , (11.33) 

d U ^ 

where D is the coefficient of self-diffusion. The solution of 
Eq. (11.33) with p(r,o) = o(t-rn) l&ads to 

-3/2 r T 

P(r,t 7 r ,0) = (STcDt) exp [- — J # 


(11.34) 
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noting the fact that v;e are here interested, not in the case of 
a molecule diffusing relative to a fixed point, but in identical 
molecules diffusing with respect to each other. The correlation 
functions of F^(t) , given in Eq. (II.25), for intermolecular 
dipolar interactions are then given by 


G^(t) 


-V- 


( 11 . 35 ) 


where n is the number of spins per molecule, M is the molecular 
weight, p is the density of the liquid, and {H and are the 
orientations of r and r respectively. Each molecule is 

/ 

considered as a hard-sphere of radius 'a' with the spin located 

-fs/ TOAct % 

at its center so that the lower limit X bhe integrals in 
Eq. ( 11 . 35 ) should be 2a instead of zero. With these considera- 
tions the correlation function becomes 


(t) = 


ny B 
(2a) 


exp(- 


u^t) 


(11.36) 


where ' 13^2 ^ Bessel function and u = 2ae where e is a 

variable entering through the Fourier transform of 
(r - r ) 2 

exp [~ — ] , The spectral density J(co) (see Eq. (11,18)) 

oDt 

is then given by 
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s rr^r _ / \ n 2 UT du ~~\ 1457 

J (qu) -p [ “I 2 2 2 -1 5 

•-o u+qo)-!-"^ ^ 

m nP rF3/2(“>]^ 

^ — ir T235-3 i ^ 


rA (^xcC) 


X 

3 


U' 


X du 

* 

4 


u- 


7 


In the extreme-narrowing limit this gives 

J(0) = 


47t n P 


(11.37) 


75 M 2aD 

If the intermolecular dipolar interaction is the only mechanism 


of relaxation Eqs, (II. 22) , (li. 2315) and (11,37) give 


r_l_] 

"^1 inter 




_ ^ Y 

5 Ha 


L. j. 


(11.38) 


The above derivation has the limitations that (i) the radial 

distribution function of the liquid has not been taken into 

account and (ii) no allowance was made for the fact that the spins 

are not necessarily at the center of the molecule. These facts 

2lB 52 

have been taken into account by several authors ' . Both the 

corre iC* tions add to the above relaxation rate and may be 

significant when taken together. These contributions could not 

be considered in a meaningful manner for the liquids studied 

due to lack of data. It is interesting to note, however, that 

the proton relaxation rate in benzene calculated with the above 

formula agrees remarkably well with the experimentally determined 
26,43 

inter* view of this we shall ignore the correction terms 

and use Eq, (II, 38) in the present analysis. 

Further, from the physical nature of the parameter ’a' 
it is clear that 'a' is not related in any simple manner to the 
molecular dimensions. It could be considered as the radius of 
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that 'equivalent sphere' of which if the liquid, is composed the 
macroscopic transport coefficients and other properties are 
reproduced. If the molecular translation in these liquids is 
viewed in terms of Stoke 's law (which is rigorously valid only 
for colloidal particles) 'a' is related to D and h through the 
relation 


kT 

&Kr]D 


(11.39) 


This reciprocal relationship between D and "H as a function of 
temperature is examined in several liquids for which the data 
on "n is available by measuring D using the spin-echo technique. 

These results/ given in Appendix indicate that the constancy <3f 

13 n btxt ^ 

■ 7 =“' obtains in several liquids, the corresponding 'a' has no direct 

X A 

relationship to molecular dimensions. 

In the present work it is necessary to have a reasonable 
estimate of 'a' for the purpose of calculating the contribution 
of intermolecular dipolar interaction to relaxation. In the 
absence of any specifically superior method it was done by using 
Eq. (11,39) whenever the data on ti is available and using the 
experimental D values. This procedure has its obvious short- 
comings in view of the discussion given above although the 
pattern of results that emerged on this basis for the molecules 
considered indicates that in these cases it has been satisfactory. 

The best method devised so far for determining intermolecular 

1 2 

contribution is by measuring H-*- and D relaxation times in 
mixtures of compounds and their deuterated analogues. This’ 
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procedure is, however, detailed and cumbersome and difficult 
to use in every case for reas o ns *^< 4 iaavai lability of compounds. 

An approximation like the present is then necessary. 

(2) Rotational diffusion; Intramolecular interactions ; 

We shall ndtr consider correlation functions involved 

in intramolecular dipolar interactions and spin-rotation 

interactions that arise from molecular reorientation on the basis 

15 

of the anisotropic rotational diffusion theory given by Huntress 
The rotational diffusion equation is given by 

|_ ^ (11.40) 

Where p(n,t) is the probability that the molecular orientation 

is given by q , , at the time t, are the principal diffusion 

constants and are the angular momentum components of the 

molecule in the principal axes system of the diffusion tensor. 

This theory exploits the formal similarity between Eq, (II,40) 

and the Schrodinger equation for a rigid rotor (with replaced 

by is the corresponding moment of inertia) . The 

conditional probability then becomes a Green's 

^ . 53 , 

function and can be written in terms of ■‘out complete ortho- 
nomal set of energy eigenfunctions of the asymmetric 

rigid rotor as 

* 

P(fio/0;a,t) = I (SJq) ' Cll.41) 

where is the eigenvalue corresponding to 
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(a) Irytramolecular dipolar interac tion -- The lattice 
functions F'^^Ct) and the spin operator for this interaction are 
given in Eqs. (11.24) to (11,27). F'^(t)’s are given in a space- 

fixed system. It is convenient to express them in a molecule- 
fixed frame for the purpose of calculating the correlation 
functions. Denoting the molecule-fixed quantities with primes 
we have 


F‘3(t) = I (-1)'*-" 1>‘2> (t)F 

-q-m 


im 


(11.42) 


m 


( 2) 

where V ^ ’s are the Wigner rotation matrices of the second 
rank. The time-dependence is now entirely contained in these 
V ^ ^ matrices and the spectral density J(o) becomes 

CO — I -.-I ^ ^ . 4r , . , . 

J(0) = J dt I (-1) 

0 ni,m' 


,( 2 ) 


(0) p'- CT,> 

-q-m <3* 


av 


(11.43) 


with 


-q-m q ' -m ' 


-q-m 


+q' -m' 


av 


(11.44) 

In Eq. (11.44) the solution of P(P,(^aO;5^ ,t:) given by Eq. (II,41) 
may be substituted. The asymmetric rotor eigenfunctions given 
in Eq, (II.41) may be expanded in terms of a complete set of 

O 1 

symmetric rotor eigenfunctions 

J ^ ^ V — J 


h,M ^ 4 

K. 


(11.45) 


where 9 p; ]yi(^5) are the symmetric rotor functions well known to 
be related to the matrices. Using the orthogonality 
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properties of the V matrices gives 

^ 2 ) / s ( 2 ) 

= A. (m)y. (-m') y. ^ ^ ^ 

c; rr, _ I a 0\ V-L-*-* 


J(o) 




£^ 2 ) 


46) 


It may be noted that the molecule-fixed frame chosen here is 
the principal axis system of the diffusion tensor. Equation 
(11.46) reduceSj^for CU dipoioM. \jA<txoJ^n, belween o- of 


J(0) 


5 


(11.47) 


Where r is the intemuclear distance in the correlation function,. 
The rotational correlation time is given for a planar 
asymmetric rotor by 


Tr 


— [ (D + D_)cos^ 9+ (Dv + D_)sin^9 — cos^«{tein^9l 

4d„ X to y s + Dc, 

R ^ . 

(11.48) 


and for a symmetric rotor by 


(6D_|_) ’^ [ 1 + 


3(D_l - D//) 
(5D^ + D//) 


_ 3(D_l- D// ) 

sin 0(1+ 

2(D^ +2D//) 


sin^e)] , 
(11.49) 


where , Dy and refer to the rotational diffusion constants 
along the x,y and z axes# D^y- and Dj^ to the diffusion constants 
for motion about and perpendicular to the symmetry axis , respect- 
ively# (e /Cp) refer to the polar angled' 'of the interriuclear' vector# 
~ ^^^x^y ^y^z ^z^x^ and Dg = ~ (D^ + + D^) . 

Substituting Eq. (11.47) in Eq, (11.22), evaluating double 
commutators and comparing with Eq, (II. 23b) we have 

1 _ 3 1i2 

^ ^l^intra DJ). 2 


1 


(11.50) 



35 


for the relaxation rate due to the dipolar interaction between 
a pair of identical spin-Jj_ nuclei. 


I t 

Spin- rotation interaction s The Fh s and s 

for the spin-rotation interaction are given in Eqs. (11.28) to 

(11.31). The spin operators are transformed to a space— fixed 

15 

frame using the Wigner p matrices of the first rank 


a'^ 

molecule -fixed 



Cl) 


( a ) 


. M. 

^space -fixed • 


(11.51) 


where Q represents the Euler angles of the molecule— fixed axes 
with respect to the space-fixed axes. Using Eqs. (11.28) to 
(II, 31) and (11.22) the equation of motion for <1^^ 
extreme-narrowing approximation can be written as 

„ m+k ^ , 

- = I ("D C(llJl,m-qqm)C(llJi,m'-q'q‘m’)C‘i^?"^ 

jl.q,k,m, ^ ^ 

q' /m' 


X / oC“C (O) ^kq^*(0) 

X j' (-m' +q' ) 5 . 51 V {[a""^, [ A^,!^]] (cT-d^)} 

(11.52) 


From Eq, (ll,52) it is clear that in the calculation of 
relaxation rate for spin-rotation interaction two types of 
correlation functions are involved viz. a rotational correlation 
function arising from the modulation of the elements of C by the 
molecular rotation and an angular velocity correlation function 
( J = I. CO where I is the moment of inertia tensor and m is the 


• ? - s 

angular velocity) which arises from the modulation of due to 
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angular velocity changes of the molecule. The separation of 
the total correlation function into a product of two implies 
the assumption that the orientation changes and the angular 
velocity changes have no significant correlation. The 
expression for T^ is now obtained to be 



SR 


m*4*Tc* rn 1 

I '-1) ■ C(llJl,m~qqm)C(llil,m'-q'q'm*)C„ 

0 

( if ( 1) 

X < 11 (o) V, , (11.53) 

Tcq k-q av 


where indicates that k takes values +1 and -1. If the angular 
velocity correlation time j is much shorter than ^5 the rotation' 
al correlation time, the significant correlation function involves 
only the angular- velocity- dependent terms, so that 

1 On i* m I 

[ ] ~ I C ( IIS' ,m-qqm) C( llJl ,m' +q-qpn) C 

”^1 SR S-,q,m,m' 

r <j'*(^--q) ( 0 ) j' (t;)> dt . ( 11 . 54 ) 

0 

Note that Eq. (11.54) is given in the principal axifes system of 
the moment of inertia tensor. It involves cross correlation 
terms between angular velocity ( w ) components. If cross- 
correlations of the type /°°<(Djji(o) <^(t) > dx with m n 

0 

(where m and n take values x,y,z) and contribution from off- 
diagonal elements of the spin-rotation interaction tensor are 
negligible and remembering that the diffusion constanf^ ^ 


D = /^„(o) ‘^(p) > dx , 

n Q “ dw 


(11.55) 
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and in the diffusional limit Eq, (11.55) gives 

('i:j)ikT 

Di r: ' 

1 


(11.56) 


where (t j) ^ is the correlation time for u and is the 
corresponding moment of inertia, the relaxation rate is given 


^ I [C^^-2Ct(Cii-C^)] Ii('Cj)i, (11.57) 

Ti SR 3 i=^,y,z 


“ T J,y,. -2Ct<=u-=t)] • 


(11.58) 


If the reorientational motion is described in terms of 

25 41 

the rotational analogue of the classical Langevin equation ' 


given by 


N = - £. 05+ n' (t) 5 

^ ^ *0 


(11.59) 


where N is the torque experienced by the molecule/ C is the 
rotational friction tensor and N' (t) is a fluctuating torque 
arising from the changes in the molecular environment, it can 
be shown that the angular velocity correlation time (neglecting 
the cross terms in the corresponding Euler's equations) is 


given by 


I. 

(Vi = -r ' 

H 


(II. 60 ) 


which relates it to the elements of the friction tensor. 


Further, it can be shown that are related to the inter- 


molecular potential through 

2Ii ^V(sN) 


^i 


(11.61) 
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where V(R^) is the potential energy of N molecules in terms 
of all the relative position vectors. Equations (II. 6o) and 
(11.61) serve to interpret the angular velocity correlation 
time. 

The Onsager criterion^^ for the motion of a molecule in 
a liquid to be diffusional is given in terms of the ratio of 
the time period for a free rotor (tf) to tj 



■'Zf 

r 





[~] 


(11.62) 


If many collisions (in which u changes) occur in one free- 


rotational cycle the reorientation is likely to be diffusional. 

It can be seen from Eqs .‘^k.dl. 5o) and (11.58) that the 
relaxation rates arising from dipolar interaction and spin- 


rotation interaction have opposite dependences on and 
therefore opposite dependences on temperature. 


II. D 


lations of Motion for Multispin Systems 


It is well known that the equations of motion of 
two spins relaxing through their mutual dipolar interaction 
leads to coupled differential equations for the z-components 
of the individual magnetizations'^' These can be 
obtained by using Eqs. (11,22) and (11.24) to (11,27) for two 
soins li and I. , to be 

- 1^^]+ -I -loj]' 

Sr r JII.63) 

a J(o)[ + 5J(o) [< -loj]. 
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where Iq 3 _ is the value of sit equilibrium. If the dipolar 

interaction is intramolecular J(0) is given by (see Eq. (II, 47)) 


J(0) = I 
5 

and if it is intermolecular in 


6 ' 
r 

origin ^'V' 


( 11 . 64 ) 


J(0) 


= 22L L 

75 Ma D 


( 11 . 65 ) 


When there are several spins in a molecule such equations 
of motion can be written for the spins involved in every intra- 
molecular and intermolecular dipolar interaction that exists. 

It should be noted that in order to write the equations for the 
intermolecular interaction it is necessary to consider at least 

^ "if 

two molecules at a time and that the correlation function is 
the same for the interaction between any two spins on different 
molecules. If the relaxation is contributed primarily by inter- 
molecular interaction the equation of motion for all the spin- 
pairs add up to give a single exponential decay for the total 


In writing these equations cross-correlation between different 
pairwise dipolar interactions may be neglected since these 
are shown to contribute a small effect^®’ . In a recent calcula- 
tion57 ^ which anisotropic effects are included, it was 
shown that these might lead to non-exponential T 3 _ -decays. From 
the experimental observation it appears that such effects are 
not significant for the present problem. 

* -k 

This implies that this interaction is operative primarily at 
the time of molecular collisions which in turn are assumed 
to be binary in character i.e. some sort of a bimolecular 
complex is formed during the collision. 
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magnetization of a multispin system. On the other hand, if the 
relaxation is only through intramolecular dipolar interactions in 
a multispin system the eq'aationsof motion for all the spin pairs 
add up to give a single exponential only if all the spins experi- 
ence intramolecular dipolar interactions of the same strength 
(e.g. protons in benzene). Normally both the types of interactions 
contribute and the relaxation decay would then be a single expo- 
nential only if all the spins experience intramolecular dipolar 
interactions of the same strength. 

For the purpose of deriving the equations of motion appro- 
priate for a given multispin system it is convenient to group 
together the spins which experience intramolecular dipolar inter- 
actions of the sSme strength. If a molecule containing N spins 
which has a such groups containing n^^ , n.2 ... .na spins so that 


n^^ + n2 + ng(= N , ( 11 . 66 ) 

the equation of motion of the total magnetization of the first 


group due to intermolecular interactions is given by 

i yfi *T' w ■ 


dt 


i=l 


j(o)< z 


i=l 





where and 1^^ represent the z-components of the spins on 
the two molecules/ J(o) is given by Eq. (11.65) . Ihe first 
term represents the interaction between n^ spins each on the 
two molecules and the other two terms represent interactions 
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between groups. These are referred to as ' cross -relaxation* 
terms. For the intramolecular interactions, however, the 
equations would be more complicated, e.g. , the equation of 
motion for the total magnetization of the first group is 


given by 


n-i 




d 

d-t 


,11.^1 - . ^11 --rl -rl 


< y I •> = y { 20 . .<1.-1 .> + ct.<i .-I. > } 

31 Z1 -^01 ^ 13 Z3 30 


i<3 


A L - "of + 


i=l ^=2 3=1 


( 11 . 68 ) 

IP 

where is the coefficient, (5/2)J(o), corresponding to the 
interaction between spin i and J on groups 1 and p respectively. 
It is clear that a sum of Egs, (11.67) and (11.68) does not 
lead to a single exponential decay for the total magnetization 

5o 

of either all the groups or any individual group. Brooks et al . 
have considered this problem and experimentally demonstrated the 
existence of two exponentials in the proton relaxation decay 
of 1-phenylpropyne. Recently similar effects were observed in 
the three proton system of 1, 1, 2-trichloroethane^® . 

In Chapters IV to VI we consider multispin systems in which 
the equations of motion of the above type yield two exponentials. 
The interpretation of the observed relaxation rates, even if it 
is q single exponential experimentally, in terms of such equations 
improves the accuracy in the determination of the relaxation 
parameters therein. 
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THE SPIN-ECHO APPARATUS * 

The phenomenon of spin-echos originally introduced by 
Hahn^^ has bean used to considerable advantage in the study 
of spin-relaxation. The response of the spin system to specially 
chosen combinations of strong resonant r. f. pulses leads to the 
formation of transient resonance signals known as 'echos'# 
characteristic of the coherence in the motion of the spins. 
Measurements on the echos for specific pulsing arrangements 
leads to the determination of relaxation times. The principal 
advantage of this method over the conventional steady-state ninr 
method is that no r.f. field is required to observe the signals 
and hence the system is completely unperturbed when the 'spin- 
echos* are observed, A brief description of the formation of 
spin-echos and details of the spin— echo apparatus and accesso- 
ries fabricated for the measurement of relaxation times and 
coefficients of self-diffusion is given in this chapter. 

Formation of Spin Echos 

Consider a system of nuclei in a magnetic field 

given by 

H = Hq? +AHQ(z)k + 2H^cos ojtr.t . (III.I) 

The first two terms represent a uniform static magnetic field 
(H^k) along the z-axis superposed by an inhomogeneity (AHQ(z)k) 

The details of the spin-echo aoparatus given in this chapter 
appear as a Technical Report 
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assumed to be along the same axis for simplicity. The last 
term represents an r.f. field of frequency close to the 
resonant value. The oscillating field is equivalent to the 
superposition of two counter-rotating fields 

^ A 

Hj^Ccosut i + sinut j) . (ill, 2) 

The precessional motion due to the first term in Eq. (III.I) 
may conveniently be described by transforming to a coordinate 
frame rotating about the z-axis at a frequency(- u k ), In this 
frame the nuclei see an effective field 

= (Hq - ^)k +AHQ(z)k + H]_ i . (III. 3) 

The effect of the other rotating field in this frame is 
negligible. If co exactly equals the resonance frequency i.e. 

= YHq , (III. 4) 

and the r.f. is sufficiently strong such that 

» aHq , (III. 5) 

the nuclei essentially precess about the x-axis. Now, if the 
r.f. field is applied just for a period t^ such that 

ti « T^ or T 2 / (III. 6) 

the relaxation terms in Eq. (l.l) may be neglected during 
this period and the nuclei precess through an angle 

e =YH^t^ , 


(III. 7) 
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about the^t-axis. Further, if before the application of the 

A 

r.f. pulse the sample attained a net magnetization 
4Q 

(Fig. lA) " this vector would also process about the "jt-axis 
through an angle 9 (Fig. IB) in the time t]_. 

Consider the case when t-|_ is such that 6 =%/ 2 . The 
entire magnetization is now tipoed into the xy-plane (Fig, 1C) 
and the r.f. is off. Then the only field seen by the nuclei 

A 

is AHQ(z)k because of which nuclei at different parts in the 

sample process at different rates. A magnetization vector 

60 

called an 'isochromat' is associated with all- the spins in 
the sample experiencing the same magnetic field. The inhomo- 
geneity results in a fanning out of these isochromats in the 
xy-plane (Fig. ID) leading to a decay of the magnetization 
in this plane in times of the order of (YAHq) This is 

known as ' f ree-induction decay' or the 'tail' of a Tc/2~pulse. 
The height of the tail is a direct measure of the magnetiza- 
tion along the z-axis prior to the application of the tc/ 2- pulse. 

The phase dispersion represented by the decay of the 
'tail' following the t/2-pulse can be reversed by applying a 
t -pulse at a time ^(^T^or T 2 ) after the t/a-pulse. The 
effect of this pulse would be to turn the isochromats through 
an angle t about a direction in the xy-plane (taken as y-axis 
for simp.licity) (Fig. IE.) and when the pulse is off the 
isochromats would continue to process in the same direction as 
they were before the % -pulse (Fig, IP) . This leads to a 
reclustering of the isochromats at a time T after the -pulse 
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giving rise to what is known as an 'echo' (Pig. IG) . The 
height of the echo is proportional to the total magnetization 
in the xy-plane at that time. 

During the time 2 t that elapsed between the initial 
■ji; /2-pulse and the appearance of the echo the relaxation terms 
in the equation of motion of the individual isochromats will 
be operative. These terms lead to a shortening of the lengths 
of the isochromats. Further loss of magnetization during this 
interval may also result from diffusion of the spins in the 
inhomogeneous field. The study of the echo for different types 
of pulsing arrangements is thus capable of yielding information 
on relaxation and diffusion. 

An apparatus to facilitate such a study was fabricated 
in our laboratory. The salient feature of this are described 
in the next section, 

iil.B Description of the Apparatus 

A spin-echo spectrometer involves an arrangement 
for applying short intense bursts of r.f. to the sample in the 
desired sequence. The induced signal after detection and 
amplification is displayed on an oscilloscope. The block 
diagram of a spin-echo spectrometer is shown in Pig. 2. The 
pulse programmer gives the required pulse sequence in the form 
of dc voltages. These gate a crystal-controlled continuous- 
wave master oscillator which gives pulsed r.f. to the coil 
surrounding the sample placed in a magnetic field. The response 
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of the sample to the disturbance is picked up by the receiver 
coil/ amplified and detected, A description of the different 
parts of the apparatus is given below. The important circuits 
are included in this chapter. Some of the routine circuits 
used are given in Appendix II. 

■k 

1. Pulse programmer ; The pulse programmer is a 
combination of Tektronix waveform generator (Tek 162) and 
Tektronix pulse generators (Tek 161 and Tek 163) together with 
some simple electronic circuits designed to give the required 
sequence of oulses. Explicit description of some of the 
important pulse sequences is given in Sec. IIIC, 

k 

2. Gated oscillator ; The gated oscillator supplies 

r.f. power of controllable amplitude and duration to the r,£, 

coil when subjected to dc pulses from the pulse programmer. 

Since at all other times|there should be no r.f, reaching the 

sample coil the carrier suppression ratio should be infinite. 

In practice this is not obtainable due to obvious limitations 

of electrical circuits. A gated oscillator with a carrier 

9 61 

suppression ratio ^ 10 designed by Blume and modified by 
62 

Proctor has been used in this apparatus (see Fig. 3. ), 

k 

We were benefitted by some consultations with Drs. K. Lalita 
and K.V.L.N. Sastry in constructing these. 

Tektronix Inc./ Beaverton, Oregon, U.S.A. 




FIG. 3: GATED OSCILLATOR 
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A crystal controlled oscillator operating at 9.78 MHz 
is used to provide continuous -wave r.f. voltage of about 

5 volts at the cathode of an isolation stage (1/2 6DJ8) . 

* * 

The crystal (.HC-6 ) is housed in an oven (HO-9) and maintained 
at 85°C by thermostat control for frequency stability. The 
isolation stage is coupled to the cathode of a loll held 
cub-off by applying (-12 v) at its plate. loll is a cylindrically 
symmetric tube whose planar grid forms a perforated region 
in the lid of a copper can which shields the entire oscillator 
circuitry. Any leakage r.f, which may reach the plate of loll 
because of its interelectrode capacitance is bypassed to ground 
through the low dynamic impedance of a conducting diode IN307, 
The dc pulses arranged in the required sequence are fed 
to one of the grids of a duo-triode (12AT7), Limitations on 
the amplitudes of the pulses obtainable from Tektronix pulse 
generators necessitate the addition of this stage to Blume's 
circuit. D.C. pulses amplified to about 90 volts swing the 
plate of loll positive. The bypass diode is back biased, 
loll conducts and r.f. of about 20 volts peak to peak reaches 
the power amplifiers through a voltage divider and an isolation 
stage (5$ 5687) . At the grid of the power amplifiers there is 
another diode which grounds any leakage r.f. Two 5763 's conneC' 
ted in parallel are used as power amplifiers to give about 300 
volts peak to peak of pulsed r.f. These tubes have a tendency 
to pick up r.f. noise and^theref ore, need good shielding, 

*From International Crystal Manufacturing Company, Oklahoma,- *• 
Oklahoma, U.S.A. 

19794 
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3» Radio-frequency receiver ; The r.f. receiver is 
an amplifier-detector system based on a circuit sent to us 
by Proctor (see Fig, 4). Pour identical stages, 6AK5's, 
each of which has an amplification factor of about lo, form 
the pre-detection amplifiers. In a spin-echo experiment one 
is usually interested in detecting signals of short duration 
6(20[isec <6 ^1 msec) and it, therefore, becomes necessary 

to broadband the receiver. Stagger-tuning the amplifier in 
pairs gives a bandwidth of about IMHz ( 3 db down) for the 
tuning range of 8-20 MHz. This places a limitation on the 
net gain of the receiver since the gain bandwidth product, 
is constant. The amplified signal is detected by a diode 
(oA5) , the output is amplified by a 12AX7 and can be displayed 
on an oscilloscope. The net gain achieved is 10,000 with a 
signal to noise ratio of 2:1 -for 80% modulation at 3|iv carrier 
voltage. Since we are using a single coil probe the receiver 
is subjected to intense r.f, pulses which saturate it. It 
is, therefore, essential to have a short recovery time if the 
actual decay is to be observed. This receiver has a recovery 
time of about 15 jAsec and accepts a maximum input signal of 
300 ^volts. 

* 

We thank Dr. W.G. Proctor, Varian Associates, Palo Alto, 

California, U.S.A. for sending us the circuits of the 

receiver and gated oscillator used by him62^ 
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To avoid oscillations the amplification stages are 
shielded from each other and the components arranged to 
minimize pick up. L 2 are high Q(=8o) coils honey-comb wound 
on a ceramic former with a ferrite core. 

The receiver was tested for linearity for output signals 
ranging from 0.5 volts to about 6.0 volts, The signal strengths 
obtained for different samples plotted against density of 
spins yield a straight line passing through the origin. 

4. Other accessories : The probe is described as 
part of the temperature variation accessory. 

'it 

The magnet used is a Varian Model V-3800 15" magnet 
with a V-FR2703 power supply. 

The signal is displayed on a Tektronix 564 Storage 
Oscilloscopel- together with Type 3A3 plug-in unit. 

The time intervals were measured using a Hewlett-Packardtih 
5245M Electronic Counter equipoed with a 5262A Time Interval 
Unit. 

"k 

We thank the Low Temperature Group of I.I.T. Kanpur for 
permitting the use of their magnet. The early testing 
of the spectrometer was done on a Polytronic (Polytronic 
Inc., Bombay, India) electromagnet loaned by Professor P.T. 
Narasimhan, Department of Chemistry, I.I.T. Kanpur. 

** 

Varian Associates, Palo Alto, California, U.S.A. 

We thank the Hydraulics Group, Department of Civil Engineering, 
I.I.T, Kanpur for loaning us this Oscilloscope. 

tt ' . . 

-Hewlett-Packard Inc., Palo Alto, California, U.S.A. 




55 


III.C Details of Performance 

1 , Typical -pulse sequences i The arrangement for 
obtaining a three pulse sequence is shown in Fig. 5. 

A Tektronix waveform generator (Tek 162) , triggered 
manually/ gives a sawtooth and a pulse (a) at the beginning 
of each sawtooth. The sawtooth is used to trigger a Tek 161 
pulse generator which in turn gives a positive pulse ('b25V) 
of adjustable width. This pulse is fed to a dif f erentiator - 
inverter circuit to obtain two positive spikes (b) and (c) . 

(a) and (b) are fed to a mixer the output of which triggers 
a Tek 163 pulse generator and (c) is fed to another Tek 163 
to obtain a set of three dc pulses. The interval between 

(b) and (c) can be adjusted by altering the width of the 

Tek 161 pulse. The interval between (a) and (bjc) can be 
varied by changing the point on the sawtooth at which Tek 161 
is triggered. The arrangement shown is for a three-pulse 
sequence in which the first two pulses have the same width 
(e.g. a sequence used for Tj_ measurement). The 

output of the two Tek 163 's are mixed and inverted and used 

to gate the master oscillator. To obtain a two-pulse sequence 
(e.g. a %/ 2-%/2 or t/S- 1 or t - tc/ 2) the initial trigger (a) is 
eliminated. 

Typical series of two and three r.f. pulses are shown 
in Figs. 6 and 7 respectively. 




fig. 5; BLOCK DIAGRAM FOR OBTAINING A THREE PULSE SEQUENCE 
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FIG. <o : A TWO PULSE SEQUENCE 

(Pulse widths: 10 fj secs and 20 jj secs; Pulse separation: 
60 fj secs. ) 



FIG. 7: A THREE PULSE SEQUENCE (Pulse widths: 10 Msec, 10 Msec, 20 Msec; 
Pulse separation: 60 Msec, 12 fjsec). 
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(2) Signals obtained ; It has been explained in 
Sec. III. A that a sequence of a pulses leads to the 

formation of an echo. Adjustments for %/2 and % pulses can 
be made using a procedure given by Schwartz" Figure 8 shows 
a tail following a t/2-pulse in glycerine. An echo, following 
a r\./2-'t-% sequence, in glycerine is shown in Pig, 9, 

III.D Measurement of Relaxation Times and Diffusion Constants 

(1) T-}_ measurement: One can employ several pulse 

sequences for measuring T-j_ e.g. Ti /2-%/2j Ti-%/2t %/2~n/2- % etc. 
The 71/2-71/2-71 method is described here. An initial 7n/2-pulse 
tips the equilibrium magnetization into the xy-plane. A 
7i;/2-i;-7c sequence applied t seconds later monitors the magnetiza- 
' tion that has developed along the z-axis during that interval. 
The height of the echo appearing a time t after the TC~pulse 
would be a measure of this magnetization and hence the time 1 
(«T;j_ or T 2 ) should be kept constant (about 0,6 msec in our 
experiments) throughout the experiment. The rate of growth of 
magnetization along the z-axis is determined by measuring the 
height of the echo as a function of t. The time constant of 
this growth gives T]_. t can be varied from a fraction of a 
second upto over a hundred seconds depending on the value of 
T 3 _ being measured. The error in T^ measurement is about + 5%, 

( 2 ) D measurement : The discussion of spin-echos in 
Sec. III. A is based on the assumption that the spins are always 
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located at toe same value of the magnetic field so that their 
precession frequencies remains constant. The height of the 
echo would then be a measure of the decay of magnetization in 
the xy-plane during the interval 2T because of relaxation. In 
actual practice, however, spins diffuse through an inhomogeneous 
field and thereby undergo changes in their processional frequen- 
cies as a function of time. There would, therefore, be an 
additional dephasing of the isochromats arising from self- 
diffusion of the molecules containing the spins thereby causing 
a decay in echo amplitudes much faster than that due to relaxa- 
tion, Carr and Purcell have shown that this decay is not an 
exponential, but has the form”^^^^ 

2-, 2 5 

M(t) = M(0)exp(- -fe- - 5IJLJL) , (111.8) 

T2 12 

in a magnetic field of constant gradient G along the z-axis, 

where M(t) is the echo amplitude at any time t. 

In all the liquids considered the exp(-t/T 2 ) decay was 

very slow because of large relaxation times and was completely 

swamped by the diffusion decay so that the first term in the 

M ( t) . 

exponential may be ignored. A plot of log 1 versus gxves 

M(0) 

D provided G is known. 

A constant gradient was obtained across the sample by 
keeping it in a homogeneous field and superposing on it a 
gradient obtained by using two circular rings of current, placed 
symmetrically with respect to the r.f. coil, with their axis 
along the homogeneous field direction. The coils were of about 
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5o turns each of No. 20 enamelled copper wire wound on perspex 

formers 'v7.o cms in diameter and placed 5,5 cms apart. 

The gradient was determined by a procedure described by 
6 5 

Douglass and McCall based on the measurement of echo shape. 
The shape of the echo at resonance for a cylindrical sample 
with its axis peirpendicular to the gradient direction is given 
by 


Jl(YGta/2) 

Y Gta/2 ' 


(III. 9) 


where 'a' is the diameter of the sample and J 2 (YGta/ 2 ) is a 
Bessel function of first order. An echo shape for a gradient 
of 0.2 gauss/cm is shown in Pig. lo. The two first minima of 
the echo shape on either side of the central maxima are separated 
by an interval t such that 

20. 54 

YG = ^ , (III. 10) 

at 


from which G may be determined. This values agrees with that 

determined from the known values of D for H 2 O®® . For the coils 

used a gradient of 1,55 gauss/cm is obtained for a current of 

loo mAs. The error in the D values measured is about + 5%. 

49 

Carr and Purcell suggested a method of determining T 2 

0^ 

by using*ji,(?t/2-pulse followed by a series of Tt -pulses. This 
largely reduces the effect of diffusion on the decay. The 
spectrometer used is constructed to make T 2 measurements also. 
Since T 2 measurements are not used in this work, the 
details of this are described in Appendix III, 
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FIG. 10 ; SHAPE OF AN ECHO IN GLYCERINE 
(Xaxis = 2 m secs/ cm. ) 
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III.E Variable Temperature Probe 

A cryostat suitable for relaxation studies from 
liquid nitrogen temperature up to 600°K was constructed on the 
basis of a design by Gutowsky^ It consists of a 
silvered glass dev/ar which contains a brass coolant chamber, 
heater coils and a probe which fit inside the dewar. A 
cross section of the assembly is shown in Pig. 11. 

The orobe consists of a split copper block 3.0 cms in 
diameter. The lower half serves merely to hold the sample. 

The upper half in addition to holding the sample supports the 
heater coils and the heat leak. The heater coils consist of 
about 2.5 meters of No. 27 double-silk-insulated Eureka wire 
wound on two cooper spools about 3 mm in diameter. These 
coils have a resistance of about 30£^ each and are connected 
in parallel. They are heated by current supplied from an 
APLAB regulated D.C. power supply. The heat leak is a 2o cm 
long brass -copper rod which forms a direct thermal contact 
between the cooper block and the coolant chamber. There is a 
10 cm long copper sleeving which can be moved over the brass 
part of the rod thereby controlling the amount of brass in 
* 

We thank Professor T.M. Srinivasan, Department of Physics, 
I.I.T. Kanpur, for suggestions regarding the construction 
of this probe. 

"k k 

Applied Electronics Ltd., Thana, Maharashtra, India 
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the circuit. Electrical connections to the power supply are 
made through a brass tubing, 5mm o.d. , which extends from the 
copper block to the top of the dewar. Another thin walled 
brass tubing vrith teflon spacers serves as an r.f. coaxial 
cable to the sample coil which lies between the two copper blocks. 
The coolant chamber is an annular brass structure with its i'nner 
and outer diameters as 3 cms and 13 cms respectively and about 
11 cms in height and can be filled with the coolant. The 
temperature is measured with a copper-cons tan ton thermocouple 
one end of which sits in the sample tube. The voltage is mea- 
sured with a Leeds and Northrup* 8687 volt potentiometer. Tempera- 
tures ranging from room temperature upwards were obtained by 
using the heater coils only. The current required for the 
maximum temperature was of the order of 1.0 amp. Re-equilibrium 
on increasing the heater current was established within 45 minutes. 
The temperature was stable to within a degree over long periods. 

III.F Sample Preparation 


The samples were contained in pyrex glass tubes 
1 cm o.d. and about lo cm in length. In all the cases the 
experiments were performed along the liquid-vapor equilibrium 
curve by filling the sample tubes in a calculated ratio of 
liquid to vapor at room temperature. This ratio is given by 


Yi 

V. 


V 


P * 

r 


(III. 11) 


Leeds and Northrup Inc, North Wal4s, Pennsylvania, U.S.A, 
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where and refer to volume of liquid and vapor at room 
temperature and and to the critical temperature and room 
temperature densities respectively. The r.f. coil essentially 
covers the region of the liquid at room temperature. A picture 
of the sample tube before sealing is shown in Fig. 12. A 
constriction was made between the liquid and vapor phases at 
room temperature to avoid the effect of diffusion on measure- 
ments. A small cavity was made near the base of the tube to 

support the thermocouple. A tube containing electrodes fixed 

★ 

with a barium getter was attached to the main body of the 
sample tube, 

'Analar' grade chemicals were first fractionally 

ic-kic 

distilled at atmospheric pressure. Dissolved oxygen was 

removed by first flushing the sample with nitrogen and then 

subjecting it to six or more freeze-pump-thaw cycles (at a 

pressure '^1 micron of Hg) and sealing at the top constriction. 

The getter was then fired following a procedure similar to 

that given by Moniz et al.^*^ The sample was then sealed at 

the second constriction. This sample was contained in the 

r.f. coil and good thermal contacts were ensured to minimize 

temperature gradients thereby reducing convection effects. 

*Obtained through the curtsey of Prof. M. Satyam of the Indian 
Institute of Science/ Bangalore and Bharat Electronic Ltd, / 
Bangalore, India. 

The chemicals thiophene and 1, 2-dichloroe thane were kindly 
provided by Professors C.N.R. Rao and M.V. George respectively 
both from the Department of Chemistry, I.I.T. Kanpur. 

'k'kk 

We thank Dr* D. Devaprabhakara and Miss Indu Mehrotra, Department 
of Chemistry, I.I.T, Kanpur for help in distilling the samples. 
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III.G Photographs of the Spectrometer 

A photograph of the spin-echo spectrometer 
console is shown in Fig. 13. Figure 14 shows a photograph 
of the magnet and cryostat assembly. 
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FIG. 13: SPIN-ECHO SPECTROMETER CONSOLE. 




fig. 14; MAGNET AND 


CRYOSTAT ASSEMBLY. 
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PROTON SPIN RELAXATION AND MOLECULAR MOTION IN THIOPHENE AND FURAN 


^"7. A Introduction 

In this chapter a study of proton spin- lattice 
relaxation and molecular motion in the liquids thiophene and 
furan is presented. The results on and D in tliiophene 
are shown in Figs- 15 and 16 respectively. Tj_ varies from 24.0 
sec at 300°K to 72.0 sec at 550°K with a maximum of 94.0 sec at 

CT CT 

490°K. D varies from 2.4 x 10 cmvsec at SOO'^K to 38.0 x lo*" 
cmvsec at 550 K. For furan (see Fig. 17) T-]_ varies from 39.0 

sec at 295 '^K to 66.0 sec at 4o5°K with a maximum of 95.0 sec at 

o — 5 2 

373 K and D varies monotonically from 5. 5 x lo"* cm /sec to 

-5 2 

36.0 X 10 cm /sec in the same temperature range. The contri- 
bution of the inte molecular dipolar interaction is separated 
from the observed relaxation rates on the basis ofja translational 
diffusion theory (see Sec. II. C) using the D values measured. The 
remaining contribution is analysed in terms of the intramolecular 
dipolar interactions and the spin-rotation interactions using 
the theory of anisotropic rotational diffusion (see Sec. II. C) 
to obtain the rotational diffusion constants of the molecule and 
the correlation times for reorientation and angular velocity 
changes and their temperature dependences. An estimate of the 
spin-rotation interaction constant of protons is made for both 
the molecules. 

* 

The work on thiophene molecule in this chapter would be 
published in J. Chem. Phys. ^ (1971)51. 

^ * "iKa AUOtk ±0 -U. X. . P 
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» ^869 * 

The geometrical arrangement * of the four protons in 

these molecules (see Fig. 18) introduces an additional feature 

into the analysis of the present experimental data which may be 

of importance in other molecules as well. This arises due to the 

fact that the two protons close to the sulphur (thiophene) and 

oxygen (furan) experience much weaker intramolecular dipolar 

interactions compared to the other two protons. It can be shown 

on the basis of Sec. II. D that the approach to equilibrium of the 

z-component of total magnetization is governed by a superposition 

of two exponentials. The decays observed in the present case, 

however, showed a single exponential (see Fig. 19) since, as the 

analysis showed, the exoonential with the longer time constant 

never has an amplitude greater than one-tenth of the other. 

IV. B Relaxation Rates in Thiophene and Furan 

The only magnetic nuclei in these molecules are the 

1 3 

four protons (apart from the 1% natural abundance of C ) . The 

occurrence of the maxima in the T^ vs lo /T graphs (Figs. 15 and 

17) shows that the spin-rotation interaction of protons is 

operative as a mechanism of relaxation in addition to the inter- 

and intra-molecular dipolar interactions. 

An important distinction appears between the relaxation 

rates of the two protons labelled I^ , I 2 and those labelled 

o 

Si , S 2 . For example, in thiophene Si and S 2 are 2.64 A apart 

o 

and are at the same distance (2,64 A) from I^ and I 2 respectively, 

o 

I^ and I 2 are 4.54 A apart. The distances are approximately the 
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same in furan. Glnce the relaxation rate arising from the 
dipolar interaction is proportional to the inverse sixth power 
of the internuclear distance the dipolar interaction between 
(Ii;l 2 )? negligible compared to the 

other three (Sj_; 32 ) and (l 2 /*S 2 ). The protonSl]_ and 

1 2 / therefore# experience much weaker intram.olecular dipolar 
interactions compared to and S 2 . 

For the purpose of writing the rate equations it is 
convenient to define 


and 



(IV. 1) 


S = S. + So . 

^ Ji ^ 


The equations of motion for < 1^> and < 82 > can be obtained, 
by using Eqs. (11.38), (11.63) to (II. 68 ) and adding the spin- 
rotation interaction, to be 




^ < I > 

dt z 






kgC <S^> 


P > + 




K 


^ * 


■[2(<S, 


> - 


(IV. 2) 


So ) + ("k" - k h. 


With the rate constants kj , kg and kg given by 

kj = (IV. 3) 

kg = Tg , 

(IV.4) 

kg = ^ , (IV. 5) 
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where 3 ^ is the distance between and S 2 ; that between 

I]_ and S]_ (or I 2 and S 2 ) ; and tjg are the correlation times 

appropriate for the dipolar interactions (S^;S 2 ) and (Ij_ 7 S^) 

(or(l 2 ;S 2 )) respectively. and are the 

relaxation rates if the spin-rotation interaction and the inter- 

molecular dipolar interactions respectively were responsible 

for relaxation. Equations (IV. 2) are identical to those given by 
50 

Brooks et al. for the relaxation of two different groups of 
nuclei on a molecule. In the present case the rate constants 
are explicitly identified in Eqs, (IV. 3) to (IV. 5). These 
authors have shown that the solution of Eqs, (IV. 2) representing 
the growth of the total magnetization =< 12+33 > following 
a t/2-pulse to equilibrium value is^ in general, a sum of 

two exponentials given by 


with 


Mo - M2 


= exp(- 

- h^t) 

+ a exp(-k t) , 

(IV. 6) 

= itn - 

kg + 

2kQ + R) / 

(IV. 7) 




(IV.8) 

= +[(kj-- 

CM 


(IV. 9) 


The relaxation decays observed in thiophene and furan, 
however, do not reveal two exponentials since a^ always has a 
value <0.1. The relaxation rate observed is then given by 
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k 


obs 


(TJ r 

1 obs 




) ] + [ (O 1 +[Y^2(i_T;^ 


-6 


s s 


IS^ 


+ -ji {-|- Y^'t^Tg r“S^ + { ^ 


(IV. lo) 


It must be noted that is not a linear superposition of the 

contributions of the type +(T^^)~^ as 

one would be tempted to write if the cross -relaxation effects 
are not examined. 


IV. C Analysis 


For the purpose of analysing the experimental results 
we shall define the observed relaxation rate as 


where 


with 


and 


^'^l^obs ^obs ^inter ^intra ' 


^intra ^SR ^DD ' 


1 r- , ^ inter. -1 

^inter ~ 5(6T^ ) 


^SR 


. SR,-1 

(T^ ) , 


(IV. 11) 

(IV. 12) 

(IV. 13) 
(IV. 14) 


= yV[ 


\)D- ’ ■■ ■- T • s ‘s 


—6 

I „ tr. + I 


r‘^1 + 

IS ^IS ^ 


inter “I ■-6\2-ih 


+ {(3Tt ) 


'’^IS ^IS 


2 2 


(IV. 15) 
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The division of relaxation rates in Eqs. (IV. 11) and (IV. 12) 
is merely for convenience and the subscripts are somewhat 
arbitrary since k^j^ter include all the contribution 

from intermolecular dipolar interactions and contains some 
cross-relaxation terms arising from, the same. 


1. Intermolecular dipolar interactions : Using the 
results of Sec, II. C we have, in the 'haird-sphere* limits the 
intermolecular relaxation rate as 


A" 


inter -1 _ 4 ^ n p 


(IV, 16) 


In the above expression p values have been calculated using 

room temperature densities and temperature coefficients for 

thiophene available in literature*^^. Viscosity data*^^ on 

thiophene between 2o°C and 80 °C used along with D values gave 

o 

a value of 1.54 + 0,05 A for ‘a' , Since the data is not available 
for furan the same value of 'a' has been used. Substituting 
the relevant constants one gets 

^inter*‘^^°P^®^®^ = 5.506 x lo"”^ (^) , 


7 p (IV. 17) 

^inter^-'^^^^^ 6.808 x 10 (^) . 

The values of / ^inter with the resulting l<^j_ntra 

for thiophene and furan are plotted in Figs. 20 and 21 respectively. 


2. Intram.olecular dipolar interactions and 
spin-rotation interactions ; 

Ihe analysis of k^^tra carried out on the 

basis of the theory of anisotropic rotational diffusion discussed 









85 


in SeCo II. C. Since both these molecules are planar and have 
an additional plane of symmetry the principal axis of the 
diffusion tensor and the moment of inertia tensor coincide^^ 

(see X and Y marked in Fig. 18) . The third axis is perpendi- 
cular to the plane of the molecule. Using Eq. (11.48) and 
the structural data of the molecules the correlation times tg 
and^Tig can be calculated. kgp, is given by Eq. (11.58). Since 
there is no information regarding the spin-rotation interaction 
tensor for different protons it will be assumed that this 
mechanism contributes equally to each of them. The values of 
or C.J- are not known for these molecules. Further analysis 
of the data on the basis of Eqs.(ii.48) and (11.58), as they 
are given, is difficult because of too many unknown quantities. 

In order to proceed with the analysis we assume that tj is 
isotropic i.e, it is the same for the three principal axes. 

This relates the three diffusion constants ini'ratios of the 
GorresDonding moments of inertia. There is no a priori justi- 
fication for this assumption. If the analysis on this basis 
yields reasonable values for the diffusion constants and their 
variation with temperature the approximation would be acceptable. 
Under this aporoximation Eq,( 11.58) becomes 
★ 

This approximation implies that the molecular collisions affect 
the reorientation about the three axes with equal efficiency. 

From Eqs. (II. 6o) and (11,61) it can be seen that this approxima- 
tion also implies that the ensemble average of the second deri- 
vative of’ the intermolecular potential is isotropicl^' '^1. 
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’^SR = ='°1 


(IV. 18) 


where 


= ' = - =t>hl . 

Which is a constant for the molecule. Substituting -Po = 90° 


and 9i 


111 3' and as 62.8568. 93.2944 and 


156.2164 amu A respectively for thiophene°° one gets 


0. 1007 


and I 


0. 1206 

D 


(IV, 20) 


For furan 9^ and 9i3 are 9o° and 165°5l‘ and 1^^ , ly and 

o2 

are respectively 54.6720, 53,5126 and 108.2304 amu A so that 


0.1153 


and ^-r 


0. 1140 


(IV, 21) 


Further, r^ and r^^ are 2.64 A for thiophene and are 2.7 56 A and 
o 

2.735 A respectively for furan. We have then 


(h ) 

DD^ thiophene 


3.301 X 10® 1.015 X 10® 

+ 


^ furan 


X [1.586 + (1 + 1.086 X lo' 

2.6763 X 10® 0.776 x lo® 

^ 

x[ 2.107 + (1 + 1.756 X 10 


-1= 


(IV. 22) 


■15 Ld ,) 2]*2 

D -L ’ 


The condition for a minimum in k^^tra considered 

in terms of Eqs. (IV. 18) and (IV. 22) by treating the quantity 
under the radical in Eq. (IV. 22) as constant with resoect to 

pD , 

temperature since the contribution of is small and does 

not vary as rapidly as D , This gives 
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Mol. 

T * 

min 


0 min 

C 'xlo^'^ 

Thiophene 

462°K 

1.30 

5.2 

3.160 

Furan 

373°K 

1. 28 

i i 

4.5 

2.852 


The values of D , and TT are reasonable. To see the implication 
of C* we shall equate it to (-^ ? ^i^ value if just 

the scalar part of the spin-rotation interaction is considered. 
This gives 


C 4 - = — I < 1.00 KHz (for thiophene) 


3 i 
43 


(IV. 23) 


l,4o KHz (for furan)* 

Powles and Figgins^'" estimate = 3.8 KHz for benzene from 
relaxation measurements. -Ramsey^^ estimates < l.o KHz for 
benzene from molecular beam data. In view of this the values 
of given in Eq. (IV. 23) seem reasonable. 

Substituting for C in the expressions for and 

using Eq. (IV. 17) allows the determination of Dj_ at any 
temperature. The values of Dj^ along with the corresponding 
Tj’s for thiophene and furan are plotted in Figs. 22 and 23 
respectively. It may be noted that the condition Tj«7:s » '’'IS 
required for writing Eq. (11.58) is satisfied. Figures 20 and 


21 also show Kgj^ and K^q as functions of temperature for the 


two cases 
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FIG. 22: TEMPERATURE DEPENDENCE OF AND FOR THIOPHE3SE. 
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IV. D Discussion 


Figure 22 shows that upto 500°K in thiophene 
obeys an Arrhenius equation: 

Dj_(T) = Dj_(o)exp [- EDj_/RT ] , (IV, 24) 

where Dj_(t) is the value of Dj^ at T°K, is the activation 

energy and R is the gas constant. The deviation at higher 
temperatures is probably due to the fact that the reorientational 
motion in the liquid is not completely diffusional and inertial 
effects come into play. The activation energy obtained from the 
plot is 3,2 kcal/mole. In the case of furan (see Fig, 23) 
does not have an Arrhenius character at all. 

For thiophene the value of XjXsee Eq. (II, 62)) at 300°K is 
about 7,0. It drops to about 2,0 at 480°K and is about 0,8 at 
SOO^K. The motion of thiophene molecules about an axis perpendi- 
cular to the molecular plane can be considered diffusional upto 
about 48o°K. For furan Xjvaries from 10.9 at 292°K to 1,1 at 
4o3°K so that one might expect the motion to be diffusional upto 
380°K. 

A computation of the amplitudes a^ of the two exponentials 
in Eq. (IV. 6) with the parameters determined by the analysis 
shows tJiat the largest value of a_ over the temperature range 
considered in both the cases is less than 0.1. Ihe longer time 
constant would, therefore, not be observable with the precision 

s. 

of our experiments. 
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The determination of the contribution of the three 
mechanisms and the parameters thereof rests on two approximations 

(a) the contribution of the intemmolecular interaction is given 
by a 'hard-sphere' model translational diffusion theory apd 

(b) the angular velocity correlation time is isotropic. 

The error introduced by (a) may not be serious since this 
approximation has been found to be valid for other liquids^^'"^^ . 
The second approximation was made to reduce the number of 
independent parameters entering the theory of rotational diffu- 
sion, The results obtained on and tj 
satisfactory. 


indicate that it is 
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PROTON SPIN RELAXATION IN M0N03R0M0THI0PHENES* 

In this chapter the experimental results on proton spin- 
lattice relaxation times in 2— and 3— bromothiophenes in the 
temperature range 300°K to 45o°K and their analysis is 
presented. The experimental results for T 2 _ and D for 2- and 
3 -bromothiophenes are shown in Figs. 24 and 25 respectively. 

V. A Analysis and Discussion 

The analysis of the relaxation data for both these 
molecules has been performed in terms of two significant 
mechanisms, namely, the inter- and intramolecular dipolar 
interactions. Spin-rotation interaction is expected to contri- 
bute very little in the temperature range studied. Measurements 
were not made at temperatures higher than 45o°K due to incipient 
decomposition. 

The molecular structures of 2- and 3 -bromothiophenes are 
shown in Fig. 26, Assuming that the interproton distances in 
these molecules are approximately the same as those in 

cr 1 

thiophene it can be seen that in 2 -bromo thiophene only the 
dipolar interactions (l];_;S) and (I 2 ;S) are significant while in 
3 -bromo thiophene only {SitS 2 ) is significant. Ihere are, in 
principle, two correlation times entering the calculation for 
2 -bromo thiophene but for the present analysis these are set 
*The work in this chapter would appear in J. Chon IPhys.^ C 197 2)’71 




s«cs 




2 5 M 

10^/T m CK)”^ 

FIG. 25: COEFFICIENT OF SELF -DIFFUSION, D, ASA FUNCTION OF 
TEMPERATURE IN 2- AND 3-^ BROMOTHIOPHENES. 






98 


the growth of the net magnetization = <1^ + > following 

a 7^/2'- pulse is a sum of two exponentials 


o z 

Mrt 


- a_j_ exp(-k^ t) + a_ exp(- k_ t) , 


(V.3) 


where = (l^ + s^) and 


\ = + R) , 


(V.4) 


with 

and 


_ ir , _ (fs - + 3^(fs - ^ 

- -[.1 + — ^ ] , 

(V.5) 

R = +[ {(kj - kg) + ak^Cfg - fj) + 4fgfj k(?]^ 

(V.6) 

The experimentally observed T^-decays in both the cases, 
however, consist of single exponentials since, as the final 
analysis showed, the longer time constant is either too long 
or has too small an amplitude to be observed. The observed 
relaxation rate is, therefore, given by 




(V.7) 


for 2-bromothiophene and 

1 . , inter. -1 •:> inter -2 , o .2 

- 5(6T, ) + 4 L4-t: + r (6T, ) + ^ 


^Ti obs 


4 ^6 


^4" 


^4 + 2 


^ Q X. X 


{V.8) 
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for 3-bromotliiophene. 

Rotational correlation times for the two molecules can 
intsiT —1 

be obtained if, (T^ )“ , given by 


, inter _i 
(T^ ) 


— 'did 

5 Ma 



(V.9) 


can be estimated at all temperatures. The D values used are 
obtained experimentally over the desired temperature range. 

A value of 'a' compatible with D and the coefficient of 
viscosity ' T) ' over this temperature range could not be 
determined since data on ' ri • is not available. It is, there- 
fore, assumed that *a* has the same value as for thiophene 
o 

viz. 1.54 A for both the molecules. This gives 


(T^nter^“l ^ 2.552 X lo“^ (§) . (V, lo) 

Density values at different temperatures were estimated using 
known densities at 2o°C and the temperature coefficients for 
thiophene. 

intBiT 

Eliminating (T^ ) estimated as above from the 

observed relaxation rates the rotational correlation times for 
both the molecules are determined over the temperature range 
studied. These are shown in Fig. 27. The two curves are in 
sharp contrast considering the fact that the molecules are 
somewhat similar. It can be seen that the variatiorj, of tj, 
with temperature for 2-bromo thiophene is similar to what is 
normally expected. The activation energy for this variation is 






lol 


2.0 kcal/mole. For 3-broinotiilophene, however, % ^ remains 
essentially unchanged with temperature indicating a large 
value for the activation energy. It must be noted that the 
values given in Fig. 27 are not very accurate since several 
assumptions were made regarding the structural and motional 
constants due to lack of data. However, the contrast in the 
trend of variation of t j- , even after allowing for the errors 
in the above estimates, is unmistakable. 
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PROTON SPIN-RELAXATION IN LIQUID 1, 2-DICHLORO ETHANE 

VI . A Introduction 

This chapter presents a study on proton spin-lattice 
relaxation in 1/ 2-dichloroethane. T^ and D have been measured 
in the temperature range 300°K to 430°K. The results are shown 
in Fig. 28. The analysis of the results on this system shows 
that the relaxation behavior and the motional characteristics 
of this molecule are somewhat different from what one might 
expect from the results on the systems described in Chapter IV 
and V. We are here, however, dealing with a molecule of 
different shape. The possibility of scalar coupling with 
chlorine nuclei and dipolar interactions modulated by ro tamer 
exchange contributing to relaxation have been examined and found 
to be negligible. The observed relaxation rates have been 
analysed in terms of inter- and intramolecular dipolar 

using the results of Chapter II and it is seen that the rotational 

/ 

correlation time Ty is essentially independent of temperature. 

It is suspected that the motion may not be diffusional and the 
possibility of inertial effects governing the molecular motion 
has been investigated. 

VI .B Analysis 

An inspection of the molecule shows that the various 
mechanisms which might contribute to relaxation are the scalar 






/ 
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7 

coupling of the protons with chlorine nuclei , intramolecular 

72 . 

dipolar interactions modulated by ro tamer exchange as well as 
rotational reorientation and the intermolecular dipolar interac- 
tions. There is no evidence for a significant contribution 
from spin-rotation interaction. Scalar coupling of the protons 
with chlorine nuclei might act as a relaxation mechanism for 

35 

the protons because of the rapid quadrupolar relaxation of Cl 
37 

and Cl nuclei. The expressions for T^ and scalar 

coupling (A I.S) of spin I with a quadrupolar nucleus S are 

7 

given by 


i ^ 2A^Si . S J^ t 12 

^ 1 +(»j. - 2 


1 

T 


2 


A^S(S+1) 

3 


{ +p } (VI,2) 

/ 2 ^ 

1 + - Ug) 


Where andr 2 'the spin-lattice and spin-spin relaxation 

times for S and uj and cog are the Larmor precession frequen- 
cies of spins I and S respectively. In the present case 
(coj - 60 X 10^ and the quadrupolar relaxation times are 

usually of the order of a microsecond or longer. For A(=2n;j) 
Where J is a few Hz in magnitude, it can be seen that the contri 
bution of this interaction to T^^ is negligible. It contributes 
only to t"^ and can, sometimes, be observed as a broadening of 
the resonance of spin I. For the discussion of T^ results this 

therefore, be considered. 


interaction need not. 



Prom the observed relaxation rates the intermolecular 

contribution can be separated by using Torrey's formula (see 

Eq. (11,38)), The density and viscosity data for this molecule 

are available in literature'^^. Using the observed values of D 

along with viscosity data gives a hard-sphere radius 'a' =1,39 + 
o 

0,05 A for this molecule. The estimated intermolecular contri- 
bution is plotted in Pig. 29. It can be seen that the tempera- 
ture variation of this interaction accounts for most of the 
temperature dependence of the observed relaxation rate. 

Among the intramolecular dipolar interactions the interac- 
tion between the protons in each of the CH 2 CI groups are the mos 
significant. (The interaction between protons on different group 
are comparatively weaker due to the larger internuclear separa- 
tion) . The molecular motion which modulates this dipolar 
interaction is likely to be rather complicated since, as is 
well known, the CH 2 CI groups undergo rotamer exchange by 
rotation about the C-C bond. The motion is, therefore, a 

hindered rotation superimposed on the molecular rotational diffu- 
73 

sion. Hubbard has considered the problem of hindered rotation 
of a CH^ group superposed on the anisotropic rotational motion 
of the molecule containing the group and calculated the relaxa- 
tion rate due to the dipolar interactions between the three 
protons. In the present case we are concerned with a single 
internuclear vector undergoing hindered rotation about the C~C 
axis. The correlation functions of angular orientation would 
include the correlation functions of the rotation of the 
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inter nuclear vector in a molecule-fixed frame and also those of 
the orientation of the molecule-fixed frame with respect to a 
laboratory-fixed frame. Such correlation functions can be 
calculated for various models for hindered rotation. Using one 
in which the CH^Cl group can assume three equally spaced angulai 
orientations about its axis of rotation and the system jumps 

"I 

from one orientation to another with a probability pet unit 

73 

time the spectral density is calculated to be 


J(0) 


5 ^6 


1 


[ 


5D^ + + 3(8ti^) 

2(D^+ 2D^) + 3(2t;j^) 



. (VI.3) 


where D (=D ) and D_ refer to the rotational diffusion constants 
along the x (in the plane of the triangle containing the two 
protons and chlorine nuclei) and z axis (C-G axis). gives the 
correlation time for rotamer exchange. 

Whatever be the nature of motion the relaxation rate from 
intramolecular interactions can be written as 


(T,) . 


1' intra 


= A X 


. 4^2 


T^eff 


where 'z 


ef f 


2 JJi ^ 

is an 'effective' correlation time, Hubbard's 


(VI, 4) 


calculation gives (from Eq. (VI. 3)) 


1 : ^°x + + 3(8 i;r) 


-1 


'eff 


6D. 


'x^2(D„ + 2D-) + 3(2To)“^ :* 


-1 1 


(VI. 5) 


X 




The correlation time for rotamer exchange can be calculated 
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using the relation 
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-1 


3a! 

2%h 


exp(- 



(VI. 6) 


where is the energy barrier for rotation. The values 

of Ea are hnown to be 5.3 and 5,1 kcal/raolc respectively for 

tbe trans and gauc he ro tamer positions. In the temperature 

range of interest ('’•'r) varies from 3.4 x lo sec at 

300 °K to about 48.0 x 10^ sec“^ at 43o°k. The observed 

-12 

values of teff are about 10 sec and hence rotamer exchange 
should give no significant contribution to relaxation. In the 
absence of rotamer exchange 


where tj. is 
given by 


( 7 


= _L 

2 -y6 


t r 


(VI. 7) 


intra 

the correlation time for reorientational motion 


1 ^^X ^'7 

’ (vi.s) 

if the motion is diffusional. From the experimental results 
on intramolecular contribution is calculated over the entire 
temperature range and is plotted in Fig. 3o. 


VI. C Discussion 


It can be seen from Fig. 30 that 'i ^ varies very 

slowly with temperature. This behaviour indicates that the 

molecular motion in this case may not be diffusional (which 

usually results in an Arrhenius type of variation with tempera- 
41 

ture) . Steele has obtained an upper limit on the value of tj^the 
correlation time for inertial type of reorientation of a 
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spherical top according to which 

r. < i T“^ (VI.9) 

^ “ 2 3k 

For the case of dichloroe thane I ^ has been replaced by 
|(l“ ^ where Ix / ly and are the principal 

moments of inertia. From known values of has been calcula- 

ted at various temperatures. These values are also shown in 
Fig. 30. These inertial correlation times are lower than the 
observed tj-'s by almost a factor of 10. However, the trend in 
the temperature dependence of Tjj^is similar to the experimental 
values . 
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Miller and Gordon have studied proton relaxation in a 
number of chloroe thanes by determining and T 2 as a function 
of concentration in OS,^ at room temperature. They also find 
that the inertial effects probably determine the rotational 
motion in these molecules. The temperature variation of the 
rotational correlation time observed in the present work also 
points to the same conclusion. It may, however, be necessary to 
study other chloroethanes as a function of temperature and also 
compare the results with a theory which is rigorously applicable 
to them before concluding definitely on the inertial effects 
in their reorientational motion. 



Chapter VII 


SUMMARY AtO CONCLUSIONS 
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SUMMARY AND CONCLUSIONS 

In this chapter some of the significant results of this 
work are summarized and discussed. 

The experiments cons is ted in the measurement of relaxation 
times as a f^inction of temperature. The analysis of these 
results involve sorting out the various mechanisms that might 
contribute to relaxation in the system and then separating the 
contribution from each mechanism. In all the cases considered 
the intermolecular interactions have been estimated on the basis 
of a theory of translational diffusion given by Torrey^"^ by 
using values of the coefficient of self-diffusion measured 
in situ using the spin-echo apparatus. This theory also involves 
a hard-sphere radius 'a' for the molecules. This radius is not 
related in any simple way to the molecular dimensions and may 
be considered as the radius of that equivalent sphere of which 
if the molecules of the liquid were composed the various 
macroscopic properties of the liquid would be reproduced. In 
the present work 'a' is estimated by using the reciprocal 
relationships between D and the coefficient of viscosity '"n.* on 
the basis of Stokes-Einstein formula. Although it does not 
seem possible to ensure that this is an accurate method the 
results indicate that in the absence of a direct method for 
determining this contribution this procedure is both useful 
and satisfactory. 



114 


Another important aspect considered, in the analysis is 
that in molecules containing groups of spins which experience 
intramolecular dipolar interactions of different strengths 
the relaxation decays may contain multiple exponentials. Further, 
the relaxation rates are not linear superpositions of contri- 
butions from the different mechanisms but are complicated 
functions involving them, because of cross relaxation effects. 
This improves the accuracy of the determination of the relaxa- 
tion parameters. 

In thiophene and furan spin-rotation interaction contri- 

butffisnsignif icantly to relaxation in addition to dipolar interac- 

tions as indicated by the maxima in vs lo /T curves. An 

analysis of the intramolecular interactions by the theory 

15 

of anisotropic rotational diffusion led to the determination 
of the rotational diffusion constants (!>,■ _) and angular 

1 —X. ^ y ^ Zj 

velocity correlation times, t ^ , as functions of temperature. 

U 

In performing the analysis r-, is assumed to be Isotropic to 

U 

simplify the expressions involved. Although there is no a priori 
justification for this assumption the results indicate that 
it is satisfactory. The rotational diffusion constant perpendi- 
cular to the plane of the molecule, , shows an Arrhenius 
type of temperature dependence in thiophene but not in furan ^ 
which is rather intriguing. It appears that in both molecules 
the reorientation departs from a diffusional character at high 

temperatures. The spin-rotation interaction constants have 

44ie 

been estimated for both^cases. The substantial amount of 
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information gathered regarding the intricacies of the molecular 
motion in these systems is indicative of the usefulness of the 
spin-relaxation method for this type of study. Furthermore/ 
this method is convenient for determining tj and estimating 
spin-rotation interaction constants of the order of a few 
kHz. 

The rotational correlation times / t , deteimained for 

r 

2— and 3-bromo thiophenes show striking contrast in their 
dependence on temperature. In the temperature range 300 ‘^K to 
45o'^K i: ^ in 2— bromothiophene shows an Arrhenius type of 
deoendence while it remains essentially constant in 3-bromo- 
thiophene. A slow variation of tj- with temperature can be 

41 

expected if inertial effects dominate molecular reorientation . 

It is, however, not clear how the motion could be so widely 

different in character for two very similar molecules. 

The temperature dependence of rotational correlation times 

in 1/ 2— dichloroethane show that inertial effects seem to govern 

the reorientation in, agreement with the conclusion reached 

72 

by Miller and Gordon . A comprehensive study in several similar 
molecules may provide somewhat conclusive evidence for the 
same. 

The spin relaxation method for the study of molecular 
motion has the limitation that it does not probe the ’shape’ of 
the -correlation functions but depends on the aisas enclosed by sijch 
functions (correlation times) , If the results obtained on the 
features of molecular reorientation in the liquids studied here 
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arean indication, it seems that this method is, nevertheless, 
capable of providing information of considerable value and 
significance in gaining insight into the liquid state. 
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APPENDIX I 


TEMPERATURE DEPENDENCE OF THE COEFFICIENT 
OF SELF-DIFFUSION IN LIQUIDS 


(Preprint of a paper by Anup 
to be published in 


Kitchlew and B.D. Nageswara Rao, 
Mol. Phys. 21 (1971).) 


» 
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TExMPERATURE DEPENDENCE OF THE COEFFICIENT OF SELF -DIFFUSION 

IN LIQUIDS 

The variation of the coefficient of self -diffusion, D, 
with temperature is of some interest in the study of the liquid 
s bate ' . In the present note measurements of D versus T are 
reported for ten different organic liquids in the temperature 
range. 300°K to 380*^K, The measurements were made using a spin- 
echo spectrometer operating at 9.78 MHz by using a J^TC - n; 

3 

seq[uence in the presence of a uniform field gradient The 

4 

details of the method and apparatus appear elsewhere , The 

error in the D values is about + 5 -per cent. The experimental 

results are shown in figures A1 and A2. 

The temperature and pressure dependence of D have been 

1-3, 5-7 

studied for several liquids by the N.M.R. method ' . Of 

the liquids chosen in the present work, to our knowledge , data 

7 8 

exists only for benzene ' . Our value of D for benzene at 

room temperature agrees well with that obtained by the N.M.R. 

8 

method and by the radioactive tracer method . The temperature 

variation reported here is, however, larger than that obtained 

o ^ 

by the latter technique in the range 15 C to 45 C, 

The Stokes-Einstein equation and also a hard-sphere 

model for liquids^ predict a simple relation between D and the 

coefficient of viscosity/T), viz. that Dt]/T is essentially 

constant. The constancy of Dh/T was verified for the case of 
1 2 

water * . For all the liquids studied values of il in the 







D X 10® cm^/ sec. 
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FIG. A2: VARIATION OF COEFFICIENT OF SELF-DIFFUSION, D, WITH 

TEMPERATURE FOR 0-, m-, p-XYLENE, NITROBENZENE AND 
THIOPHENE. 
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temperature range of interest are available in the literature^^. 
The quantity Dri/T at various temperatures was calculated using 
interpolated values of D and T], These are given in the table. 
For all but the last two liquids this quantity remains fairly 
constant. The activation energy for translational diffusion/ 

Ed = - R ["Sjjn D/s(l/T)] , calculated using figures^l andA2/ is 
also given in the table. In cases where Jln D versus 1/T plots 
are not straight lines Eq is an average value. 

If the average value of D'H/T is equated to k/6Tta(k is the 
Boltzmann constant) the hard-sphere radius 'a' compatible with 
the Stokes-Einstein relation may be obtained. These radii for 
the first eight liquids are given in the table. For the three 
xylenes / benzene and chlorobenzene/ the ’a' values are nearly 
the same. In all cases 'a' is smaller than the molecular dimen- 
sions and for pyridine and piperidine it is much smaller. If we 
use the expression given by Longuet-Higgins and Pople Dti/t = 
Skc/lOTia' where c is the fractional volume occupied by the 
molecules in the liquid and a' is the hard-sphere radius related 
to 'a' used above, through a' = (9c/5) (a) . The value of c is not 
known in general, but even for c 1, a' is still smaller than 
the molecular dimensions in most cases. While the variation of 
D and 11 seems to follow the temperature dependence predicted by 
the Stokes-Einstein equation and the hard-sphere model in most 
cases, the apparent hard-sphere radii seem not to have physical 
significance in either theory. 





table 


. No. 

Liquid 



Dr]/T X 

10^° 


a 

o 

(A) 

kcal/mo] 

27°C 

40°C 

60°C 

o 

o 

o 

GO 

ioo°c 

1. 

Benzene 

5.47 

5.95 

6.45 

6.32 

6. 22 

1. 20 

3.5 

2. 

Ch lorobenz ene 

6.40 

7.15 

8.26 

8.32 

3.10 

0.96 

3.7 

3. 

Pyridine 

25.8 

23.0 

21.0 

20. 1 

19.1 

0. 34 

2.1 

4. 

Piperidine 

32. 1 

25.6 

22.3 

22.0 

23.0 

0. 29 

1.6 

5. 

Thiophene 

4.98 

4. 12 

4.66 

4« 66 

4.83 

1, 58 

2.7 

6. 

p- Xylene 

6.35 

5.76 

7.01 

7.13 

6. 69 

1.08 

3.3 

7. 

o- Xylene 

5.80 

5.74 

5.70 

5.92 

6.34 

1.24 

3.2 ' 

8. 

m- Xylene 

5.40 

5.40 

5.52 

5.88 

6.38 

1. 28 

3.2 

9. 

Ni t robenz ene 

19.4 

17.0 

14.7 

12.6 

11.6 

- 

1.9 

0. 

Dime thy lani line 

1.75 

2.02 

2.52 

3.30 

4.43 

- 

6.3 


128 


REFERENCES 


1. J.H. Simpson and H.Y. Carr, Phys . Rev., Ill , 120l (1958) 

2. G.B. Benedek and E.M. Purcell, J. Chem. Phys. 7 ^, 2003 (1954) 

3. D..’W. McCall, and, DuG.' Douglass, J. Phys. Chem., _6_2, 1102 (19J 

4. A. Kitchlew and 3.D. Nageswara Rao, Technical Report No, 1/70, 

Department of Physics, Indian Institute of Technology, 
Kanpur (U.P.) , India(l97o) 

5. D.C. Douglass, D.-W. McCall and E.W. Anderson, J. Chem. Phys., 

34, 152 (1961) 

6. D.W. McCall, D.C. Douglass and E.W, Anderson, Phys, Fluids, 

2 , 87; (1959); Ber. Bunsenges, 67 , 336 (1963) 

7. D.W. McCall, D.C. Douglass and E.W, Anderson, J. Chem. Phys., 

1555 (1959) 

8. P.A. Johnson and A.L, Babb, Chem. Rev, 387 (1956) 

9. H.C. Longuet-Higgins and J.A. Pople, J. Chem. Phys., 25, 884 

(1956) 

10. International Critical Tables, Vol. 2 (McGraw-Hill Bobk Co) 
(1926) 



APPENDIX II 


ADDITIONAL DETAILS OF THE SPIN-ECHO APPARATUS 
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FIG./^: LAYOUT OF THE OSCILLATOR 



FIG.4^: layout OF THE GATE 
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FIG.ftS': LAYOUT OF THE RECEIVER 




FIGM: power supply FOR THE GATED OSCILLATOR AND PULSE PROGRAMMER 
(excluding the Tektronix units) 
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FIG.A^i: PULSE INVERTER 



APPENDIX III 


MEASUREMENT OF T2 BY CjVRR -PURCELL METHOD 
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MEASUREMENT OF T 2 BY CARR-PURCELL METHOD 

A standard way of measuring T 2 would be to disturb the 
equilibrium magnetization by a ii/2-pulse which tips the entire 
magnetization into the xy-plane. An: -pulse applied at an 
interval t after the ^l^/2 pulse gives an echo (at a time 2t) 
the height of which would be a measure of the magnetization 
in the xy-plane at that time. Varying t would indicate the 
rate of decay of this magnetization. 

It has already been explained in Sec. III.D that since 
the spins are located in an inhomogeneous field the relaxation 
decay would be superimposed by an additional dephasing of the 
isochromats because of diffusion. The decay following a 
^/2-t sequence in the presence of a constant field gradient, 

G, is of the form^^^ (see Eq, (III, 8)) 

4- Dv^G^t^ 

M(t) =M(o)exp(- — 1 ), (1) 

Tg 12 

Carr and Purcell^ have suggested a method of obtaining an echo 
envelope which reveals the natural decay due to relaxation even 
in the presence of diffusion. In contrast with the method 
described earlier, in this method the entire decay is obtained 
in a one-shot experiment, A series of tc - pulses at interval 
2 t from each other, are applied at a time t after an initial 
'Ji;/2-pulse. At the times 2 t, 4x.., etc., there will be echos 
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because the effect of each it -pulse would be to reverse the 

direction of the dephasing isochromats . Carr and Purcell have 

shown that the effect of interpolating a large number (n) of 

% -pulses in an interval of duration t is to reduce the mean 

square phase dispersion by a factor 1/n^- Thus the decay is 

1 2 

given by an expression of the type * 

M^p(t) = M(0)exp ^ .) , (2) 

and at any instant t = 2nx 

M(^p(t) =M(o)exp( - ~ - I— , 9 . ^ ( 3 ) 

T2 3 

A plot of M^pCt) versus t would, therefore, be an exponential 
inspite of diffusion. For small t’s the time constant of the 
decay is essentially T 2 . 

Experimental Arrangement for Obtaining a Carr-Purcell Sequence 

An experimental arrangement for obtaining a Carr -Purcell 
sequence, due to Meiboom is shown in Fig. A, 12. A switch S 
normally closed is pressed to give a voltage step which 
triggers a Teh 162 waveform generator which in turn gives 

(a) a single pulse at the moment the voltage step is applied, 

(b) a recurrent sawtooth as long as the switch is kept pressed 

The single pulse obtained when the switch is pressed is 
used to trigger the Tc/2-pulse generator and the recurrent 
sawtooth triggers the u-pulse generator. The two kinds of 
pulses are mixed and used to gate the r.f, oscillator. The 
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interval between the U /2 and the first it—pulse is adjusted 
to be half that between the subsequent n: -pulses by triggering 
the Tc -pulse generator at the midpoint of the sawtooth. A 
typical Carr-Purcell decay for proton resonance ia glycerine is 
shown in Fig, A, 13, 
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FIG. ^13; CARR-PURCELL SPIN-ECHOS IN GLYCERINE 
(t= 0.25 msecs; T 2 = 21.5±2.0msecs.) 



